AD-784  940 

EXTENDED  ARRAY  EVALUATION  PROGRAM. 
SPECIAL  REPORT  NUMBER  12.  CONTINUED 
EVALUATION  OF  THE  NORWEGIAN  LONG- 
PERIOD  SEISMIC  ARRAY  -  FINAL  REPORT 

Philip  R.  Laun,  et  al 

Texas  Instruments,  Incorporated 


Prepared  for: 

Advanced  Research  Projects  Agency 
Air  Force  Technical  Applications  Center 


3  December  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 

ALEX(OI)  STR  73  12 


EVALUATION  OF  THE  NORWEGIAN  LONG-PERIOD  SEISMIC  ARRAY  -  FINAL  REPORT 


SPECIAL  REPORT  NO.  12 

EX  TENDED  ARRAY  EVALUATION  PROGRAM 

Prepared  by 

Philip  R.  Laun,  Wen- Wu  Shen,  William  H.  Swindell 

TEXAS  INSTRUMENTS  INCORPORATED 
Equipment  Group 
Post  Office  Box  6015 
Dallas,  Texas  75222 


Prepared  for 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
AFTAC  Project  No.  VE  LA  T/ 2705/ B/  ASD 
Alexandria,  Virginia  223 14 

Sponsored  by 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Program  Code  No.  2F10 
ARPA  Order  No.  1/14 


3  December  1  973 

Acknowledgment:  This  research  was  supported  by  the  Advanced 

Research  Projects  Agency,  Nuclear  Monitoring  Research  Office 
under  Project  VELA-UNIFORM,  ard  accomplished  under  the  tech¬ 
nical  direction  of  the  Air  Force  technical  Applications  Center 
under  Contract  No.  F33657 -72  -C-0725. 

'  i 

II 


Equipment  Group 


ABSTRACT 


This  report  describes  the  final  results  of  the  continued  evalua¬ 
tion  of  the  Norwegian  Long-Period  Seismic  Array  (NORSAR)  by  Texas  Instru 
ments  Incorporated  at  the  Seismic  Data  Analysis  Center  over  the  period 
1  April  1973  to  30  September  1973. 

The  major  areas  of  study  were: 

•  Noise  analysis 

•  Array  processing  performance 

•  Effectiveness  of  matched  filters 

•  Surface  wave  detection  capability 

•  Performance  of  standard  surface  wave  discriminants 

A  total  of  133  Eurasian  events  and  36  noise  samples  were 
processed  and  analyzed  during  this  period.  The  results  were  combined 
jvith  earlier  data,  when  applicable,  in  order  to  maximize  the  data  base. 

Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  necessarily  representing  the  official  policies,  either  expressed 
or  implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Tech 
nical  Applications  Center,  or  the  US  Government. 


TABLE  OF  CONTENTS 


SECTION 


I. 


II. 


III. 


IV. 


V. 


TIT  LE 


ABSTRACT 
INTRODUCTION 
DATA  BASE 

NOISE  ANALYSIS 

A.  INTRODUCTION 

B.  DATA  BASE 

C.  TIME  VARIABILITY  OF  THE 
AMBIENT  NOISE  LEVELS 

D.  NOISE  DIRECTIONALITY 

E.  SPATIAL  COHERENCE 

F.  CONCLUDING  REMARK 

ARRAY  PROCESSING  PERFORMANCE 

A.  INTRODUCTION 

B.  DATA  BASE 

C.  MCF  AND  BEAMSTEER  NOISE 
REJECTION 

D.  SIGNAL  DEGRADATION  OF  MCF 
AND  BEAMSTEER  PROCESSORS 

E.  MCF  PROCESSING  OF  SIGNALS 
USING  ONE  NOISE  SAMPLE  FOR 
MCF  DESIGN 

MATCHED  FILTER  PERFORMANCE 

A.  INTRODUCTION 

B.  PROCESSING 


iv 


PAGE 

iii 

1-1 

II-  1 

III-  1 
III-  1 
III- 1 

III  -  2 
III -6 
III- 8 

III-  11 

IV-  1 
IV-1 
IV- 1 

IV -2 

IV  -5 


IV-  10 

V- l 
V-l 
V-l 


IM 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPU  TING.  FORM 

1.  REPORT  number 

2  GOVT  ACCESSION  NO. 

3  RECIPIENT'S  CATALOG  NUMBER 

/ID- Wo 

*  TITLE  fa.irt  Subtltlr) 

EVALUATION  OF  THE  NORWEGIAN  LONG- 
PERIOD  SEISMIC  ARRAY  -  FINAL  REPORT 

S  TYPE  of  REPORT  4  PERIOD  COVERFO 

Special 

6  PERFORMING  ORG  REPORT  NUMBER 

ALEX(01)-STR-7  3-12 

7.  AuTHORfJj 

Philip  R.  Laun,  Wen-Wu  Shen,  and 

William  H.  Swindell 

8  CONTRACT  OR  GRANT  NUMBERfM 

F33657-72-C-0725 

9  PERFORMING  ORGANIZATION  name  and  aooress 

Texas  Instruments  Incorporated 

Equipment  Group 

Dallas,  Texas  7  522T, 

10  PROGRAM  ELEMENT.  PROJECT  TASK 
AREA  »  WORl  UNIT  NUMBERS 

VELA  T/  2705/  B/  ASD 

11  CONTROLLING  OFFICE  NAME  ANO  ADORESS 

Advanced  Research  Projects  Agency 

Nuclear  Monitoring  Research  Office 

Arlington,  Virginia  22209 

12  REPORT  DATE 

3  December  1973 

13  NUMBER  OF  PAGES 

134 

U  MONITORING  AGENCY  NAME  A  AOORESSf//  dllletent  Irom  Controlling  Ollier) 

Air  Force  Technical  Applications  Center 

VELA  Seismological  Center 

Alexandria,  Virginia  22314 

IS.  SECURITY  CLASS  (ol  this  report) 

UNC  LASSIFJED 

15  t  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16  distribution  statement  (oi  ihi *  Report) 

APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 

17.  DISTRIBUTION  STATEMENT  (ol  the  ebilrect  entered  In  Bloch  30.  II  dlllerent  Irom  Report) 

18  SUPPLEMENTARY  NOTES 

ARPA  Order  No.  1714 

19  KEY  WORDS  (Continue  on  reverse  side  II  necessary  and  Identity  hy  block  number) 

Norwegian  Seismic  Array  Seismic  event  detection 

(NORSAR)  Evaluation 
Seismic  noise  analysis 
Chirp  matched  filler 


Seismic  event  discrimination 


JO  ABSTRACT  (Continue  on  reverse  side  II  nece s&ery  and  identily  hv  block  number) 

This  report  describes  the  results  of  the  continued  evaluation  :  f  the 
Norwegian  Long-Period  Seismic  Array  (NORSAR)  by  Texas  Instruments 
Incorporated  at  the  Seismic  Data  Analysis  Center  over  the  period  1  April 
1 97  3  to  30  Septembe  r  1 97  3.  "'"mfl’wAL  TECHMICAL 

The  major  areas  of  study  were: 

*  Sptinp.ftpld  VA  ?.?151 

• _ Noise  Analysis _ I 


UNCLASSIFIED 


SECURITY  CL  ASSIFIC  AT  ION  of  THIS  p  AC,  L  rut  .•>,  Pnf«  I  ntm 


DO  ,  1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


_ UNCLASSIFIED _ 

security  CLASSIFICATION  or  THIS  PAGt(m>an  Data  Fnlarad) 


20.  continued 

•  Array  Processing  Perforn  ice 

•  Effectiveness  of  Matched  Filters 

<#  Surface  wave  detection  Capability 

•  Performance  of  Standard  surface  wave  discriminants. 

A  total  of  133  Eurasian  events  and  36  noise  samples  were  pro¬ 
cessed  and  analyzed  during  this  period.  The  results  were  combined 
with  earlier  data,  when  applicable,  in  order  to  maximize  the  data 
base. 


)  '  CO 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEftFhen  Data  Fnlrrrd) 


f 


TABLE  OF  CONTENTS 
(continued) 


SECTION 


TITLE 


PAGE 


REGIONALIZATION  BY  CHIRP 

FILTER  LENGTH 

CHIRP  FILTER  SNNR  IMPROVE¬ 


MENT 


l 

.  + 


VIII. 


NORSAR  LONG -PERIOD  SURFACE  WAVE 
DETECTION  CAPABILITY 


BEHAVIOR  OF  STANDARD  DISCRI¬ 
MINANTS 

A.  M  -  mL  MEASUREMENTS 

s  b 

B.  AL-AR  VERSUS  m. 


CONCLUSIONS 

A.  MAJOR  RESULTS 

B.  SUGGESTIONS  FOR  FUTURE 
ANALYSIS  OF  NORSAR 


VI-1 


VII-1 
VII- 1 
VII-5 


VIII- 1 
VIII- 1 


VIII- 5 


REFERENCES 


IX-1 


APPENDIX  A 


APPENDIX  B 


APPENDIX  C 


APPENDIX  D 


LIST  OF  FIGURES 


FIGURE  TITLE  PAGE 

1-1  SITE  DIAGRAM  OF  THE  NORSAR  LONG- 

PERIOD  ARRAY  1-4 

H-1  m  DISTRIBUTION  OF  EVENTS  ANALYZED 

A't*  NORSAR  LCNG-PERIOD  ARRAY  APRIL 
1971  -  AUGUST  1973  H-2 

HI- 1  HISTOGRAM  OF  SITE  AVAILABILITY  III-3 

HI -2  THE  FREQUENCY  OF  LOWER  MICRO- 

SEISMIC  PEAK  OF  NOISE  SAMPLES  III-4 

III_ 3  AVERAGE  SITE  RMS  LEVEL  OF  NOISE 

(20-40  SECONDS)  HI-5 

III-4  NOISE  SOURCE  AZIMUTHS  IH-7 

HI- 5  MULTICHANNEL  COHERENCE  OF  THE 

7  AUGUST  1972  NOISE  SAMPLE  III -9 

HI- 6  MULTICHANNEL  COHERENCE  OF  THE 

4  NOVEMBER  1972  NOISE  SAMPLE  III- 10 

III- 7  MULTICHANNEL  COHERENCE  OF  THE 

8  JANUARY  1973  NOISE  SAMPLE  III- 12 

IV.  1  FULL  ARRAY  MCF  IMPROVEMENT  IN 

NOISE  REJECTION  AND  BEAMSTEER  OUT¬ 
PUT  LEVEL  VERSUS  NOISE  SAMPLE  (OFF- 
DESIGN  NOISE  0.  025  -  0.  059  Hz  or  17-40  sec)  IV-6 

IV - 2  REDUCED  ARRAY  MCF  IMPROVEMENT  IN 

NOISE  REJECTION  AND  BEAMSTEER  OUTPUT 
LEVEL  VERSUS  NOISE  SAMPLE  (OFF-DESIGN 
NOISE  0.  025  -  0.  059  Hz  or  17-40  sec)  IV-7 

IV _3  ^TiO  OF  BEAMSTEER  TO  MCF  NOISE  OUT- 

p\^  (0.  025  -  0.  059  Hz)  IV  -  1  3 

IV ,4  THbs,  RELATIVE  MCF/ BS  SIGNAL- TO- 

NOIS*  RATIOS  (0.  025  -  0.059  Hz)  IV- 15 

V- l  OPTIMUM  CHIRP  FILTER  LENGTHS  (IN 

SECOND^)  FOR  LOVE  WAVES  AT  NORSAR 

FOR  THE  ^ASSBAND  0.  025  -  0.  059  Hz  V-3 


FIGURE 
V  -2 


V- 3 

VI- 1 

VI-2 

VI-3 

VI-4 

VI-5 

VI- 6 

VII-  1 

VII -2 

VII- 3 

VII-4 

D-l 


LIST  OF  FIGURES 
(continued) 


TITLE  1  PAGE 

■  » 

t 

OPTIMUM  CHIRP  FILTER  LENGTHS  (IN 

SECONDS)  FOR  RAYLEIGH.  WAVES  AT 

NORSAR  FOR  THE  PASSBa!nD  0.  025  - 

0.  059  Hz  V -4 

SEISMIC  REGIONS  DERIVED  FROM  CHIRP 

FILTER  ANALYSIS  V-5 

NORSAR  LP  SURFACE  WAVE  DETECTION 

STATISTICS  FOR  ALL  EARTHQUAKES 

1971  -  1972  VI-3 

NORSAR  LP  SURFACE  WAVE  DETECTION 

STATISTICS  FOR  THE  KURILE-KAMCHATKA 

AREA  VI-4 

NORSAR  LP  SURFACE  WAVE  DETECTION 
STATISTICS  FOR  THE  CENTRAL  ASIA  AREA  VI-5 

NORSAR  LP  SURFACE  WAVE  DETECTION 
STATISTICS  FOR  THE  PRESUMED  EXPLO¬ 
SIONS  VI-7 

NORSAR  LP  SURFACE  WAVE  DETECTION 
STATISTICS  FOR  THE  SUMMER  EARTH¬ 
QUAKES  VI-8 

NORSAR  LP  SURFACE  WAVE  DETECTION 
STATISTICS  FOR  WINTER  EARTHQUAKES  VI-9 

RAYLEIGH  M  -  m  PLOT  FOR  EVENTS 

APRIL  1971  -PN0VF:MBER  1972  VII-2 

LOVE  M  -  m  PLOT  FOR  EVENTS 

APRIL  1  §7 1  -  NOVEMBER  1972  VII-3 

AR  -  m  DISCRIMINANT  COMBINED  1971 

AND  1972  EVENTS  VII-6 

AL  -  m  DISCRIMINANT  COMBINED  1971 

AND  1972  EVENTS  VII-7 

LOCATION  OF  IT  A  EVENTS  PROM  PDE 
INFORMATION  D-3 


Vll 


LIST  OF  FIGURES 
(continued) 


FIGURE 

TITLE 

PAGE 

D-2 

RAYLEIGH  WAVE  GROUP  VELOCITIES 

FROM  ITALY  MEASURED  AT  NORSAR 

D-5 

D-3 

LOVE  WAVE  GROUP  VELOCITIES 

FROM  ITALY  MEASURED  AT  NORSAR 

D-6 

D-4 

SPECTRAL  RATIOS  NORMALIZED  TO 

25  SECONDS 

D-7 

D-5 

GROUP  VELOCITY  CURVES  FOR  ITALIAN 
EVENTS  AT  NORSAR 

D-9 

D-6 

GROUP  VELOCITY  CURVES  FOR  ITALIAN 
EVENTS  AT  NORSAR 

D-10 

viii 


LIST  OF  TABLES 


TABLE  TITLE  PAGE 

IV- 1  IMPROVEMENT  IN  NOISE  REJECTION  BY 

MCF  PROCESSOR  FOR  THE  FULL  AND 
RF  jUCED  NORSAR  ARRAY  IV-3 

IV -2  SIGNAL  AMPLITUDE  DEGRADATION  OF 

BEAMSTEER  PROCESSOR  FOR  NORSAR 
FULL  ARRAY  AND  ITS  REDUCED  ARRAY  IV-9 

IV-3  SIGNAL  AMPLITUDE  DEGRADATION  OF 

RAYLEIGH  WAVE  FROM  MCF  AND  BEAM- 

STEER  PROCESSOR  FOR  NORSAR  FULL 

ARRAY  AND  ITS  REDUCED  ARRAY  IV-ll 

IV - 4  MCF  PROCESSED  EVENT  LIST  IV-12 

V- l  CHIRP  MATCHED  FILTER  IMPROVEMENT 

IN  SIGNAL+NOISE/ NOISE  RATIO  V-7 

A- 1  LIST  OF  EVENTS  USED  FOR  EVALUATION  A-3 

C-l  1971  NORSAR  NOISE  SAMPLES  C-2 

C-2  1972  NORSAR  NOISE  SAMPLES  C-3 

C-3  1973  NORSAR  NOISE  SAMPLES  C-5 

D-l  EVENT  PARAMETERS  AND  OPTIMUM 

CHIRP  LENGTHS  D-2 


SECTION  I 


s 

I 

I 

I 

I 


1! 

a 

i 

i 

i 

i 

i 


INTRODUCTION 


This  report  presents  the  cumulative  and  final  results  of  a  study 
of  the  Long  Period  Norwegian  Seismic  Array  (NORSAR).  The  study  began  on 
April  1971  for  the  purpose  of  evaluating: 

•  The  array  detection  capability  for  Eurasian  events 

•  The  performance  of  various  discriminants  at  NORSAR  for 

Eurasian  ’’vents 

•  Methods  cf  sustaining  or  enhancing  these  capabilities. 

These  three  objectives  were  achieved  by  the  following  studies: 

•  Noise  analysis 

•  Signal  analysis 

•  Array  processing  effectiveness 

•  Matched  filtering  performance 

•  Detection  threshold  estimation 

•  Behavior  of  standard  discriminants. 

The  study  consists  of  three  reports:  Texas  instruments  Special 
Report  No.  5,  1972;  Texas  Instruments  Special  Report  No.  7,  1973;  and  (his  re¬ 
port  which  is  the  third  and  final  report. 

The  noise  analysis  is  covered  in  some  detail  in  both  Special 
Report  No.  5  and  No.  7,  however,  some  additional  data  are  presented  here. 

The  signal  analysis  for  the  full  array  was  completed  in  Special 
Report  No.  5  and  results  for  the  partial  array  are  presented  here. 
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The  array  processing  performance  study  is  continued  through¬ 
out  the  three  reports.  The  first  report  was  concerned  with  the  feasibility  of 
routine  MCF  processing,  noise  stationarity ,  anti  signal  degradation  by  the 
beamsteer  processor.  The  second  report,  with  additional  data,  continued 
the  comparison  of  MCF  and  beamsteer  performance  for  both  the  full  array  and 
for  various  subarrays,  investigated  the  azimuthal  dependence  of  array  gains,  and 
compared  signal-to-noise  ratio  (SNR)  improvements  of  the  MCF  and  beamsteer 
processors  for  several  seismic  events.  This  report  contains  additional  partial 

array  MCF  and  beamsteer  processing  and  a  further  investigation  of  the  ef¬ 
fectiveness  of  MCF'e  designed  from  noise  which  occurs  many  days  before  or 
after  the  time  of  application. 

The  matched  filtering  performance,  detection  threshold  es¬ 
timation,  and  the  behavior  of  standard  discriminants  studies,  which  depend 
upon  a  large  ensemble  of  events  for  reliable  estimates  have  been  reported  in 
preliminary  fashion  in  the  first  two  reports.  In  this  report,  these  earlier  data 
have  been  combined  with  33  additional  events  from  August  and  November  1972 

and  the  final  results  are  presented. 

The  data  used  in  this  report  are  discussed  in  Section  II. 

Sections  III  through  VII  contain  the  various  studies  listed  above.  Section  VIII 
summarizes  the  results,  presents  conclusions  and  lists  areas  for  further  study. 
Appendix  A  contains  the  list  of  all  events  used  in  the  evaluation,  Appendix  B 
discusses  the  events  which  have  been  dropped  from  the  data  base,  Appendix  C 
contains  a  list  of  the  noise  samples  used,  and  Appendix  D  contains  some  obser¬ 
vations  on  linear  group  velocity  ch’rp  filter  design. 

The  long  period  NORSAR  is  an  array  of  22  seismometer  sites 

spread  over  a  circular  area  approximately  100  km  in  diameter,  and  is  located 
north  of  Oslo,  Norway.  Each  site  consists  of  three  orthogonal  seismometers 
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(vertical,  north-south,  east-west)  with  a  response  centered  at  25  seconds.  A 
diagram  of  the  array  is  shown  in  Figure  1-1. 

Long-period  NORSAR  data  is  received  at  the  Seismic  Data 
Analysis  Center  (SDAC)  from  Norway  by  a  communication  system  called  the 
Trans-Atlantic  Link  (TAL).  These  data  are  recorded  on  magnetic  tape,  along 
with  ALPA  and  LASA  long-period  data,  and  saved  for  future  analysis.  The 
analysis  was  performed  at  the  ‘'DAC  using  an  off-line  array  evaluation  software 
package  developed  by  Texas  Instruments  under  Contract  No.  F33657-69-C-1063. 
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FIGURE  1-1 

SHE  DIAGRAM  OF  THE  NORSAR  LONG-PERIOD  ARRAY 
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SECTION  II 
DATA  BASE 


The  dita  used  in  this  report  consist  of  the  signals  and  the 
noise  sample?  analyzed  in  the  two  previous  reports,  Special  Reports  No.  5 
and  No.  7,  plus  an  additional  36  noise  samples  from  about  two  to  seven  hours 
long  recorded  from  August  1972  to  March  1973,  133  events  from  August  and 
November  1972,  and  all  1973  presumed  explosions  through  August  1973. 

Appendix  A  contains  the  total  list  of  events  processed,  and 

Appendix  C  contains  the  list  of  noise  samples  edited  along  with  the  various 

paramete  s  associated  with  them.  A  histogram  of  the  bodywave  magnitude 

(m  )  distribution  of  the  events  processed  is  given  in  Figure  II- 1. 
b 

With  a  few  exceptions,  only  those  events  originating  in  Europe 
and  Asia  were  selected.  Since  one  of  the  primary  objectives  of  the  analysis 
was  to  obtain  good  estimates  of  detection  threshold,  many  events  were  needed 
at  the  lower  magnitudes  with  4.  5.  The  >e  events  were  obtained  mostly 

from  the  LASA  and  NORSAR  bulletins  with  the  remainder  obtained  from  the 
PDE  bulletin.  The  various  event  information  sources  are  detailed  in  Appendix 


The  events  in  Appendix  A  are  515  events  successfully  processed 
from  a  list  of  839  events  considered.  The  events  not  processed  were  rejected 
for  the  following  reasons: 

Lack  of  data  or  gaps  on  tape  175 

Parity  errors,  spikes,  other  tape  problems  17 

Mislocated  events  4 


Interfering  events 

Events  later  shown  to  be  deep  (>50  km) 


HQ- 


SJU0A3  jo  aoqumtvj 


ARRAY  APRIL  1971-AUGUST  1973 


[ 


Thus,  aside  from  the  technical  difficulties  of  recording  which 
account  for  21%  of  the  total  events,  one  can  expect  a  certain  loss  of  signals 
due  to  interfering  events  (17%  in  this  case). 

The  basic  signal  processing  consisted  of  editing  the  event 
including  at  least  1200  seconds  of  noise  preceding  the  P-wave  arrival, 
rotating  the  components  to  vertical,  transverse  and  radial,  and  beam¬ 
steering  in  the  source  direction  using  3.  5  km/ sec  for  the  vertical  and 
radial  and  4.0  km/sec  for  the  transverse.  Processing  of  the  1971  events 
was  done  using  three  different  passbands  to  determine  which  was  the  most 
suitable.  On  the  basis  of  this  study,  a  17  -  40  second  band  (0.025-0.059  Hz) 
was  decided  upon  and  used  for  the  1972  and  1?73  data.  This  particular 
band  was  not  used  for  the  1971  events,  so  1971  results  from  the  17  -  50 
second  (0.020-0.059  Hz)  band  were  used. 
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SECTION  III 
NOISE  ANALYSIS 


A.  INTRODUCTION 

The  noise  analysis  presented  in  this  section  is  an  extension  of 
the  noise  field  study  performed  in  the  preceding  years  (Special  Report  No.  5, 
1972;  Special  Report  No.  7,  1973).  The  objective  of  this  analysis  is  to  char¬ 
acterize  the  noise  field  at  NORSAR  to  improve  signal  enhancement  and  detec¬ 
tion  performance. 

The  long-period  noise  investigation  includes: 

•  Time  variability 

•  Noise  directionality 

•  Spatial  coherence. 

B.  DATA  BASE 

In  addition  to  the  noise  samples  of  Special  Report  No.  5  and 
No.  7,  a  total  of  36  new  noise  samples  were  edited  for  this  report  covering 
the  period  from  August  1972  to  March  1973.  The  duration  of  these  samples 
ranged  from  about  two  to  seven  hours  with  an  average  length  of  about  four 
hours.  A  list  of  the  complete  noise  ensemble  used  !n  this  study  is  in  Ap¬ 
pendix  C. 

The  data  were  sampled  at  two-second  intervals  and  were  re¬ 
corded  in  256-second  segments.  After  quality  checking  the  data,  cross-power 
matrices  were  then  computed  at  64  frequencies  from  0.  0  to  0.  25  Hz  for  each 
noise  sample.  Quality  checks  included  plots  of  the  vertical  components  of 
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three  or  four  sites  to  check  for  unreported  surface  waves,  and  inspection 
of  the  individual  segment  powers  and  the  site  auto -power  spectra  for  unusual 
values. 

The  average  number  of  sites  available  was  nineteen,  an  increase 
of  four  from  last  year.  Figure  III- 1  shows  a  histogram  of  the  number  of 
times  a  site  was  deleted  versus  the  site  number. 

C.  TIME  VARIABILITY  OF  THE  AMBIENT  NOISE  LEVELS 

The  comparison  of  earlier  noise  spectra  with  the  spectra  of  the 
latest  data  showed  that  the  vertical  and  horizontal  components  of  the  noise 
field  have  continued  to  show  the  same  spectral  behavior.  In  addition  the 
horizontals  contained  both  Love  and  Rayleigh  wave  energy  and  were  slightly 
less  coherent  than  the  verticals  (Special  Report  No.  5),  so  for  these  reasons, 
the  study  was  done  only  of  the  vertical  component. 

In  general,  the  noise  spectra  showed  two  miuroseismic  peaks: 
a  lower  one  at  about  0.  06  Hz  (16  second  period)  and  an  upper  one  at  about 
0.12  Hz  (6-8  second  period).  The  former  occurs  near  or  within  the  signal 
processing  band  (0.025  to  0.059  Hz  or  17  to  40  seconds)  while  the  latter 
is  far  beyond  it.  For  the  purpose  of  signal  processing  and  detection,  the 
noise  in  the  signal  band  is  our  primary  concern.  Figure  III- 2  shows  the 
frequency  of  the  lower  micros  ’ismic  peak  for  each  sample.  The  peak  fre¬ 
quency  ranged  from  0.051  to  0.083  Hz  (about  1  2  -  20  seconds)  with  per¬ 
haps  a  slight  tendency  to  be  nearer  the  lower  frequency  range  in  the  winter 
than  in  the  summer.  This  loweiing  of  the  peak  frequency  to  the  edge  of 
the  signal  processing  band  during  winter,  coupled  with  an  increase  in  noise 
level  and  coherence  (discussed  below)  suggests  that  the  use  of  MCF  processing 
may  be  advantageous  during  the  winter  months.  This  is  discussed  further  in 
Section  IV. 

Monitoring  the  seasonal  changes  in  noise  levels  was  continued 
from  Special  Report  No.  7.  Figure  III- 3  shows  the  average  single  site  RMS  level 
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HISTOGRAM  OF  SITE  AVAILABILITY 
(36  NOISE  SAMPLES) 


(20  -  40  second  period  band)  of  each  noise  sample  versus  time.  This  figure  in¬ 
cludes  the  entire  noise  ensemble  from  May  1971  to  March  1973.  These  data 
confirm  the  seasonal  noise  trends  of  the  20  to  40  second  period  band  at  NORSAR 
as  discussed  in  Special  Report  No.  7.  From  early  spring,  about  February  or 
March,  NORSAR  has  low  (7rryz  RMS) , mildly  fluctuating  noise  levels  until  late 
fall.  Beginning  in  October,  the  noise  starts  to  rise  rapidly  and  erratically  to 
the  mid-winter  months  of  December  to  February,  during  which  the  noise  level 
is  highly  variable  and  probably  dependent  on  the  severity  of  the  North  Atlantic 
weather.  In  February  or  March  the  noise  level  rapidly  assumes  its  summer 
Tevel  when  the  severe  winter  storms  apparently  abate. 

In  sunmary,  the  conclusions  of  this  study  are: 

•  The  frequency  of  the  lower  microseismic  peak  lies  in  the 
range  of  0.05-0.08  Hz  (20  -  12  seconds)  with  a  strong  ten¬ 
dency  of  dropping  to  0.06  Hz  (17  seconds)  in  winter. 

•  The  20  to  40  second  RMS  noise  level  was  about  6  to  7  m/z  in 
the  summer,  about  14  rryz  in  late  autumn  and  early  winter,  with 
erratic  fluctuations  in  the  late  winter  up  to  a  maximum.  The 
maximum  observed  level  of  79  rry z  would  increase  the  detection 
threshold  from  summer  to  winter  about  1.  2  Ms  units. 

D.  NOISE  DIRECTIONALITY 

The  directionality  of  the  noise  was  investigated  by  means  of  high 
resolution  f-k  spectra  computed  for  the  vertical  component  of  each  noise 
sample  at  the  frequency  of  the  lower  microseismic  peak. 

Figure  III  -  4  shows  primary  noise  directions  of  each  sample. 

The  solid  dots  indicate  the  azimuths  of  maximum  power  density  at  a  velocity  of 
3.  5  km/sec  for  the  corresponding  frequencies  shown  in  Figure  III  -  2.  The  range 
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of  azimuths  where  the  power  was  within  6  dB  of  the  peak  power  is  shown  by 
vertical  straight  lines  for  each  noise  sample.  Figure  III-4  shows  that  the 
noise  power  was  generally  isotropic  in  the  summer,  but  was  primarily  from 
the  north  to  northwest  in  the  winter  with  azimuths  near  0  ,  250  ,  and  280 
occurring  most  often.  The  noise  directions  of  the  last  four  noise  samples  (in 
March  1973)  were  typical  of  winter  noise  even  though  the  noise  levels  had  drop¬ 
ped  back  to  summer  values.  This  shows  more  evidence  of  the  low  correlation 
between  noise  level  and  direction  observed  previously. 

E.  SPATIAL  COHERENCE 

Multichannel  coherences  were  computed  for  the  7  August,  4 
November,  1972  and  8  January  1973  noi?  e  samples  (Figures  III-5  ,  III-6  and 
HI  -  7) .  Because  of  ‘hr*  inferest  in  small  array  (  <8  sites)  performance,  two 
cases  of  multichannel  coherence  were  computed:  first,  site  1  was  predicted 
from  all  the  other  sites  in  the  array,  and  second,  site  1  was  predicted  from 
only  the  inner-ring  sites. 

For  all  three  samples,  the  coherences  of  both  the  inner  ring 
and  the  full  array  are  nearly  the  same  although  the  inner  ring  coherence  is 
slightly  less  due  to  the  fewer  number  of  sites  used.  These  figures  also  illus¬ 
trate  some  of  the  problems  involved  in  MCF  design.  Coherence  plots  in  Special 
Report  No.  5  of  summer  and  winter  data  indicated  that  the  highest  coherences 
were  at  the  lower  microseismic  peak,  however,  Figure  III- 5  shows  that  the  noise 
from  7  August  1972  has  its  highest  coherence  at  0.08  Hz  when  the  microseismic 
peak  from  Figure  III  -  2  is  at  0.055  Hz.  At  this  latt-r  frequency,  the  coherence 
is  very  low  so  that  MCF  processing  would  not  be  worthwhile.  Figure  III-6,  for 
4  November  1972,  shows  high  coherences  within  the  signal  band  but  that  the  peak 
power  is  at  0.066  Hz,  outside  the  signal  band.  These  are  indications  of  potentially 
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FIGURE  III- 5 

,TICH  AN  INTEL,  COHERENCE  OF  THE  7  AUGUST  1972  NOISE  SAMPLE 


MULTICHANNEL  COHERENCE  OF  THE  8  JANUARY  1973  NOISE  SAMPLE 


high  MCF  effectiveness.  The  8  January  1973  sample,  Figure  III  -  7 ,  has  both 
microseismic  peak  power  and  h;gh  coherencies  outside  the  signal  band.  Even 
though  this  noise  sample  is  typical  winter  noise,  the  low  signal  band  coherence 
indioates  that  an  MCF  would  probaoly  be  ineffective.  As  will  be  shown  in  Section 
IV,  the  predictions  of  MCF  effectiveness  for  these  noise  samples  are  borne  out. 

F.  CONCLUDING  REMARK 

The  thirty-six  long  duration  noise  samples,  covering  the  periods 
from  summer  1972  through  early  spring  1973,  have  been  found  to  have  time  va¬ 
riability,  and  directionality  consistent  with  previous  results,  therefore,  their 

major  characteristics  can  be  considered  well-defined.  The  coherences  of  the 

\ 

noise  is  not  as  well  defined  ,  but  the  coincidences  of  the  microseismic  peak 
power  and  maximum  coherence  in  the  winter  gives  strong  evidence  of  discrete 
noise  generating  sources.  MCF  effectiveness  becomes  a  function  of  the  position 
of  the  microseismic  peak  in  relation  to  the  signal  processing  band.  In  the  summer 
the  peak  power  often  occurs  in  the  signal  band  but  the  noise  is  largely  incoherent 
and  isotropic  so  that  MCF  processing  may  not  be  justified. 
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SECTION  IV 


ARRAY  PROCESSING  PERFORMANCE 

A.  INTRODUCTION 

The  principal  objectives  in  this  extension  of  the  NORSAR  array 
processing  performance  evaluation  were: 

•  Continued  exploration  of  the  feasibility  of  routinely  design¬ 
ing  and  applying  multichannel  filters  (MCF's). 

•  Comparison  of  the  effective  noise  rejection  and  signal  de¬ 
gradation  of  the  full  NORSAR  array  and  a  reduced  array  con¬ 
sisting  of  only  the  inner  ring  sites. 

These  objectives  were  accomplished  by: 

•  Est’mating  the  noise  rejection  achieved  by  the  MCF  and  beam- 
steer  (BS)  processors  for  the  full  array  and  the  reduced  array. 

•  Measuring  the  signal  degradation  caused  by  the  MCF  and  BS 
processors  for  both  arrays. 

•  MCF  processing  of  signals  using  only  one  noise  sample  for 
MCF  design. 

B.  DATA  BASE 

Twenty-  three  samples  of  noise  data  for  the  period  from  7  August 
1972  to  20  February  197  3  were  used  for  MCF  design.  Both  the  MCF  and  beam- 
steer  processor  were  then  applied  to  two  portions  of  a  noise  sample  called  on- 


design  noise  and  off-design  noise.  The  former  was  in  the  MCF  design  gate 
and  the  latter  outside  of  the  design  gate.  Each  portion  was  4096  seconds  long. 

Twelve  events  from  the  1972  summer  event  ensemble  were  used 
for  signal  beamforming  loss  and  noise  reduction  measurements  for  the  full 
NORSAR  array  and  the  reduced  array.  Similar  measurements  were  made 
on  seven  August  1972  events  processed  by  an  MCF  beamformer. 

Using  the  19,  200  -  second  data  from  the  7  August  1972  noise 
sample,  MCF's  were  designed  and  applied  to  the  noise  preceding  31  events 
occurring  between  3  August  and  18  August  1972.  The  purpose  was  to  compare 
the  noise  reduction  achieved  by  MCF  and  BS  processors  on  data  increasingly 
displaced  in  time  from  the  MCF  design  interval.  Relative  signal-to-noise 
ratios  were  computed  on  the  19  events  having  moderate  signal-to-noise  ratios 
on  the  single  site. 

C.  MCF  AND  BEAMSTEER  NOISE  REJECTION. 

The  MCF's  were  computed  for  only  the  vertical  component  of 
the  data,  using  the  design  parameters  of  Special  Report  No.  7: 

•  Dispersive  signal  model  oriented  to  a  beam  direction  of  90° 

•  Signal-to-noise  ratio  equal  to  four  at  all  frequencies 

•  Two  percent  white  noise  added  to  the  data 

•  Frequency  domain  design  by  transforming,  cross -multiplying 

and  stacking  256-second  segments  at  all  frequencies  0.00- 

0.  246  Hz.  The  data  were  not  hanned 

On  the  average,  20  sites  were  available  for  the  full  array  and 
seven  sites  for  the  reduced  array. 
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Table  IV-1  presents  the  relative  MCF/BS  noise  rejection  for 
the  full  and  reduced  arrays.  Among  the  three  passbands  computed,  the  best 
improvement  in  noise  rejection  by  the  MCF  generally  was  obtained  in  the  wider 
band  of  0.02-0.  10  Hz.  For  the  reduced  array,  the  MCF  was  4  dB  more  effective 
than  the  beamsteer  processor  in  the  signal  processing  band,  however,  this  ad¬ 
vantage  dropped  to  less  than  2  dB  with  the  full  array.  Although  the  array  gain 
of  the  full  array  is  4.  5  dB  (10  log  (number  of  sites)  )  greater  than  the  small  array, 
the  use  of  an  MCF  with  the  smaller  array  could  reduce  this  site  difference  to 
less  than  1  dB. 

Previous  analysis  of  only  winter  noise  using  the  full  array  did 
not  show  much  correlation  between  beam  or  single  site  noise  level  and  improve¬ 
ment  in  array  gain.  With  additional  data  covering  the  transition  from  summer 
to  winter  noise  conditions,  the  winter  rise  in  noise  level  is  accompanied  by  a 
rise  in  MCF  improvement  although  sample  to  sample  correlation  of  level  and  im¬ 
provement  is  6 till  poor  (Figure  IV-1).  However,  for  the  reduced  array,  there  is 
a  strong  correlation  between  the  BS  level  and  the  relative  MCF/BS  noise  rejection 
from  sample  to  sample  (Figure  IV-2),  particularly  in  the  signal  processing  band. 

Multichannel  filter  array  gain  improvement  showed  distinct  changes 
related  to  the  seasonal  noise  changes.  For  the  full  array,  MCF-over-BS  gains 
in  the  signal  band  fluctuated  around  the  zero  dB  level  before  the  12  November 
noise  sample,  but  showed  significant  improvements  ra  ging  from  0  to  7  dB  with 
an  average  of  about  4  dB  in  the  winter.  This  implies  that  a  decrease  of  0.  2  M 

s 

units  in  the  detection  threshold  could  be  achieved  at  NORSAR  using  an  MCF  pro¬ 
cessor  in  the  winter.  Even  better  gains  were  obtained  for  the  reduced  array; 
these  ranged  from  0  to  9  dB  for  both  summer  and  winter  at  NORSAR  with  slightly 
less  seasonal  change  being  observed. 
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D.  SIGNAL  DEGRADATION  OF  MCF  AND  BEAMSTEER  PROCESSORS 

The  signal  amplitude  degradation  of  the  beamsteer  processor 
was  measured  using  twelve  1972  winter  ensemble  events  and  sixteen  1972 
summer  ensemble  events.  Seven  of  the  summer  events  were  also  MCF  pro¬ 
cessed  for  comparison.  Both  the  full  array  and  the  reduced  array  were  used 
for  these  measurements. 

Events  were  selected  io  have  at  least  a  moderate  signal-to-noise 
ratio  on  both  the  reference  site  (usually  site  1)  and  the  output  beam  traces.  Am¬ 
plitude  degradation  of  the  unfiltered  signal  was  measured  by  taking  the  ratio  of 
the  amplitude  of  a  particular  wave  cycle  from  the  reference  site  to  the  amplitude 
of  the  same  cycle  of  the  output  beam.  Three  such  ratios  were  averaged  for  each 
trace  and  the  average  ratio  was  then  expressed  in  decibels. 

Table  IV-2  shows  the  signal  loss  of  the  beamsteer  processor  for 
the  full  and  reduced  arrays.  Negative  values  indicate  that  the  output  beam  am¬ 
plitudes  exceeded  the  reference  site  amplitude.  Ra/leigh  wave  signal  loss  ranged 
from  -0.  8  to  6.  2  dB  for  the  full  array  and  -0.  6  to  3.  1  dB  for  reduced  array  with 
averages  of  2.0  dB  and  1 . 4  dB  respectively.  The  beamforming  loss  of  the  Rayleigh 
and  Love  waves  are  not  significantly  different.  The  smaller  loss  of  the  reduced 
array  is  significant,  however,  and  reflects  the  high  signal  similarity  among  the 
inner  ring  sites.  Loss  contributed  by  errors  in  the  estimated  phase  velocity  used 
for  beamforming  is  also  less  for  the  reduced  array. 

The  average  signal  degradation  obtained  for  these  events  are 
approximately  1.0  to  1.  4  dB  larger  than  those  obtained  from  the  1971  events. 

The  reasons  for  this  are  not  clear.  The  measurement  technique  of  the  earlier  re¬ 
sults  was  slightly  different,  using  only  one  wave  cycle  instead  of  three.  This  has 
contributed  a  part  of  the  difference.  A  second  possibility  is  that  one  of  the  channels 
had  inverted  polarity  causing  partial  cancellation  of  the  signal  for  a  few  of  the  events. 
However,  periodic  checks  of  the  data  were  made  to  check  for  that  condition. 
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SIGNAL  AMPLITUDE  DEGRADATION  OF  BEAMSTEER  PROCESSOR 
FOR  NORSAR  FULL  ARRAY  AND  ITS  REDUCED  ARRAY 


Event 

Name 

Signal  Degradation  (dB) 

Full  Array 

Reduced  Array 

No.  of 
Sites 

Rayleigh 

Wave 

Love 

Wave 

No.  of 
Sites 

Rayieigh 

Wave 

Love 

Wave 

SIN*002*10 

19 

0.  3 

0.  5 

7 

1.8 

1.0 

TAI*004*12 

17 

0.  3 

0.  7 

6 

0.  3 

0.4 

TAI*006*06 

16 

2.  5 

1.4 

7 

1.6 

0.4 

RA*006*09 

16 

2.  5 

1.6 

7 

1.8 

1.4 

SIB*013*17 

19 

1.2 

2.  9 

6 

0.  8 

2.  0 

SIN*042*05 

18 

(-0.8) 

1.  9 

8 

(-0.6) 

1.9 

SIN*  047*2  3 

16 

0.7 

- 

7 

0.  6 

“ 

KUR*057*02 

18 

2.  0 

1.  3 

8 

1.2 

0.  6 

YUN*057*18 

17 

1.  9 

1.4 

6 

2.  1 

2.  1 

TIB*075*06 

20 

2.7 

4.2 

8 

2.  0 

4.  9 

KUR*077*07 

19 

3.2 

3.5 

8 

0.  8 

1.4 

TAD*077*09 

19 

1.6 

3.  9 

8 

0.  2 

3.  8 

WRS*2  04*05 

18 

2.  1 

-1.  7 

7 

2.  1 

-0.  9 

TIB*204*1 6 

17 

0.7 

2.  1 

7 

0.7 

0.  3 

KUR*209*00 

16 

3.9 

1.2 

7 

1.9 

0.  9 

KUR*2 11*21 

19 

3.  9 

0.  8 

6 

1.  5 

0.  6 

KAM*2I6*I2 

18 

1.7 

3.  2 

8 

1.  1 

1.8 

IRA*2 16*22 

18 

1.7 

- 

7 

1.  o 

- 

TUR*2 1 6*2 1 

19 

0.  9 

1.  5 

8 

0.  9 

1.5 

TUR*217*05 

19 

2.  0 

1.  6 

8 

2.  0 

1.4 

IRA*22I*I9 

14 

4.  0 

5.  6 

6 

3.  0 

4.2 

BAI*222*1 9 

12 

0.  3 

1.  0 

7 

0.  7 

-0.2 

ERS*222*20 

17 

6.2 

3.  2 

6 

2.4 

- 

MON*23 1*  12 

13 

2.  1 

2.  0 

7 

0.  5 

- 

KUR*231*21 

14 

3.  2 

- 

7 

3.  1 

CHI*243*1 5 

18 

0.  5 

1.  3 

7 

0.  8 

0.  5 

SIN*243* 17 

18 

2.4 

6.  4 

7 

2.4 

5.  3 

CHI*243*18 

18 

0.  6 

2.  1 

7 

1.  0 

0.  6 

Average 

2.  0 

2.  2 

'  1.4 

1.6 

IV -9 


Table  IV-3  gives  the  vertical  Rayleigh  wave  signal  loss  from 
the  MCF  processor.  The  beamsteer  losses  are  included  from  Table  IV-2  for 
comparison.  Average  MCF  signal  losses  for  the  full  and  partial  arrays  were 
3.  2  dB  and  2.  3  dB  respectively.  This  is  close  to  the  average  beamsteer  loss  for 
these  events  of  2.9  dB  and  1.8  dB  respectively. 

The  event  ERS*222*20  had  a  very  high  Rayleigh  wave  loss  of  6.  2 
dB  for  the  full  array  but  a  near-average  loss  of  2.  4  dB  for  the  partial  array. 

This  event  was  found  to  have  significant  off-azimuth  (probably  multipath)  energy. 

The  Love  waves  from  this  event  showed  no  unusual  behavior. 

E.  MCF  PROCESSING  OF  SIGNALS  USING  ONE  NOISE  SAMPLE  FOR  MCF 
DESIGN 

Earlier  results  (Texas  Instruments,  Special  Report  No.  7,  1973) 
have  shown  that  a  single  noise  sample  might  be  useful  for  design  of  multichannel 
filters  to  be  applied  to  events,  days  or  weeks  before  or  after  the  noise  sample  it¬ 
self.  This  implies  that  either  the  noise  field  is  generally  stationary  or  that  it 
changes  slowly,  perhaps  gradually  reducing  the  effectiveness  of  the  MCF.  This 
possibility  would  permit  routine  MCF  processing  of  the  data  without  the  consi¬ 
derable  expense  in  computer  time  needed  to  estimate  cross-power  spectral  ma¬ 
trices.  To  further  investigate  the  possibilities  of  this  idea,  31  events  occuring 
between  the  third  and  the  eighteenth  of  August  1972  were  processed  using  MCF's  de¬ 
signed  from  the  7  August  1972  noise  sample.  These  events  are  listed  in  Table 
IV-4.  The  MCF  design  noise  sample  had  typical  summertime  characteristics  of 
low  level  and  isotropy. 

Figure  IV-3  shows  the  MCF/BS  relative  noise  reduction  in  the 
0.  025-0.  0 59Hz  signal  band  computed  in  a  1500-second  gate  prior  to  the  P-wave 
arrival  of  each  event.  Among  the  thirty-one  events,  sixteen  favored  the  MCF 
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TABLE  IV- 3 

SIGNAL  AMPLITUDE  DEGRADATION  OF  RAYLEIGH  WAVE 
FROM  MCF  AND  BEAMSTEER  PROCESSOR  FOR  NORSAR 
FULL  ARRAY  AND  ITS  REDUCED  ARRAY 


Event 

Name 

Signal  Degradation  (dB) 

Full  Array 

Reduced  Array 

No.  of 
Sites 

MCF 

BS 

No.  of 
Sites 

MCF 

BS 

TUR*2 16*21 

19 

0.  8 

0.  9 

8 

1.4 

0.  9 

IRA  *2 16  *22 

18 

2.4 

1.  7 

7 

1.  6 

1.0 

TUR*217*0  5 

19 

1.2 

2.  0 

8 

1.6 

2.0 

IRA*221*19 

14 

2.  8 

4.  2 

6 

2.  8 

3.0 

ERS*222*20 

17 

6.  6 

6.  2 

6 

3.0 

2.  4 

MON*23i  *1 2 

13 

4.0 

2.  1 

7 

1.7 

0.  5 

KUR*231*21 

14 

4.7 

3.  2 

7 

4.  1 

3.  1 

Average 

3.  2 

2.  9 

2.  3 
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RATIO  OF  BEAMSTEER  TO  MCF  NOISE  OUTPUT 
(0.  025  -  0.  05  0  Hz) 


and  fifteen  the  beamsteer.  Although  not  shown,  other  passbands  also  were  com¬ 
puted.  In  general,  the  results  are  in  agreement  with  those  shown  in  Table  IV-1. 
That  is,  the  MCF  was  superior  in  wider  bands  with  higher  frequencies,  but  not 
in  narrow,  lower  frequency  bands.  Figure  IV- 3  shows  that  the  MCF/  BS  re¬ 
lative  gains  are  rather  randomly  distributed  positively  and  negatively  and  that 
they  do  not  seem  to  be  a  function  of  the  event-noise  sample  time  separation  in¬ 
dicating  that  the  supposed  necessity  of  using  a  noise  sample  close  to  the  signal 
may  not  always  be  valid. 

Nineteen  of  these  events  had  good  signal- to -noise  ratio  (SNR) 
and  were  suitable  for  comparison  of  MCF/BS  SNR  enhancement.  Figure  IV-4 
shows  the  relative  MCF/BS  SNR’s  for  the  0.025-0.059  Hz  band.  Eecause  signal 
loss  for  the  MCF  and  BS  are  about  the  same,  the  SNR's  were  mostly  affected  by 
noise  reduction  and  consequently,  the  results  were  close  to  that  of  the  noise 
analysis  shown  in  Figure  IV- 3. 


The  results  of  this  investigation  show  that  although  the  summer 

noise  is  well-behaved,  its  almost  isotropic  character  allows  an  MCF  little  ad- 

0 

vantage  over  a  conventional  beamsteer  processor.  On  the  other  hand,  the  directional 
nature  of  winter  noise  which  offers  definite  MCF  improvements,  is  complicated 
by  severe  level  changes  and  other  transient  phenomena  which  would  effectively 
limit  the  use  of  fixed  MCF's.  Time-adaptive  multichannel  filters  would  seem 
to  offer  strong  potential  for  eliminating  the  directional  noise  while  coping  with 
the  transient  level  changes  and  should  be  investigated  in  any  future  analysis. 
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THE  RELATIVE  MCF/BS  SIGNAL- TO-NOISE  RATIOS 


SECTION  V 

MATCHED  FILTER  PERFORMANCE 


A.  INTRODUCTION 

The  goals  of  the  matched  filter  performance  study  were  to  de¬ 
termine  the  available  signal-plus-noise-to-noise  ratio  (SNNR)  improvements 
of  the  matched  filters  as  compared  to  a  bandpass  filter  and  to  regionalize  the 
events  by  matched  filter  parameters. 

The  matched  filters  used  in  previous  studies  have  been  linear 
group  velocity  chirp  filters  (linear  with  frequency)  and  reference  waveform 
filters  (RVVF).  Because  RWF's  did  not  show  a  clear-cut  superiority  over 
chirp  filters  in  SNNR  gain  (Special  Report  No.  5,  1972)  and  there  was  no  a 
priori  means  of  selecting  "good"  reference  waveforms  as  filters  (Special 
Report  No.  5,  1973)  this  report  is  solely  concerned  with  linear  group  velocity 
chirps.  Appendix  D  contains  a  study  of  the  limitations  and  errors  involved 
using  this  linear  frequency  design. 

B.  PROCESSING 

A  suite  of  five  chirp  filters  and  one  bandpass  filter  (0.  020-0.  059 
Hz  or  17-50  seconds  in  1971;  0.  025-0.  059  Hz  or  17-40  seconds  in  1972  and  1973) 
were  routinely  applied  to  the  three  component  beamsteered  traces  of  each  event. 
The  lengths  of  the  five  chirps  were  determined,  at  first  empirically  and  later 
from  the  maps  of  optimum  chirp  lengths  described  in  Special  Report  No.  7. 

The  SNNR  of  a  signal  was  measured  as  the  ratio  of  zero-to- 
peak  signal  amplitude  to  the  RMS  noise  amplitude.  The  signal  was  picked  as  the 
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maximum  amplitude  within  a  signal  gate  determined  by  a  velocity  range  of  4.  0 
km/ sec  to  2.  5  km/ sec  and  the  noise  was  measured  over  a  1200  second  gate 
usually  immediately  preceding  the  P-wave  arrival.  The  optimum  chirp  length 
was  the  length  which  produced  the  highest  zero-to-peak  amplitude. 

C.  REGIONALIZATION  BY  CHIRP  FILTER  LENGTH 

An  additional  96  LR  vertical  and  93  LQ  transverse  optimum  chirp 
filter  lengths  were  added  to  the  data  of  Special  Report  No.  7  and  were  used  to 
update  the  maps  of  optimum  LR  vertical  and  LQ  transverse  chirp  lengths.  Se¬ 
veral  events  which  were  listed  in  Special  Report  No.  7  as  having  non-optimum 
chirp  lengths  were  re -analyzed.  Most  of  these  events  were  from  Taiwan,  with 
the  rest  from  central  Asia  and  the  Mediterranean  area. 

Figures  V-l  and  V-2  show  the  results  of  the  entire  ensemble 
contoured  on  the  Eurasian  Continent.  During  contouring  it  was  assumed  that 
the  chirp  length  would  increase  monotonically  with  distance  from  NORSAR.  The 
data  also  were  smoothed  to  show  regional  variations  rather  than  point  to  point 
variations.  The  contours  have  changed  little  from  Special  Report  No.  7.  The 
region  of  the  Caspian  Sea  still  shows  the  extraordinarily  long  chirp  lengths  which 
cannot  be  reconciled  with  the  monotonicity  rule. 

The  regionalization  is  suggested  by  changes  in  the  direction 
and  gradient  of  the  contours.  Figure  V-3  shows  vhe  seismic  regions  derived 
from  the  chirp  length  maps.  These  regions  are  essentially  unchanged  from 
Special  Report  No.  7.  As  Appendix  D  shows,  the  optimum  chirp  length  is  not 
only  a  function  of  dispersion  but  also  of  the  signal  spectrum  received  at  the  array. 
The  signal  spectrum  usually  has  been  modified  by  the  effects  of  radiation  pattern, 
multipath  interference,  and  attenuation.  While  these  effects  are  important  con¬ 
siderations  in  any  regionalization  scheme,  they  have  been  neglected  here.  Signal 
spectrum  effects  are  most  likely  the  cause  of  the  unusual  chirp  behavior  around 
the  Caspian  Sea. 
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D.  CHIRP  FILTER  SNNR  IMPROVEMENT 

The  average  chirp  filter  SNNR  improvement  was  computed 
for  the  events  in  each  of  the  seismic  regions  shown  in  Figure  V-3.  The  indi¬ 
vidual  chirp  gains  for  each  event  within  each  region  are  given  in  Table  V-l. 
Included  in  Table  V-l  are  the  best  chirp  lengths,  SNNR  of  ihe  bandpassed  beam, 
the  improvement  of  the  chirp  filter  SNNR  over  the  bandpassed  beam  SNNR,  and 
the  averages  of  chirp  filter  improvements  by  region.  Non-detection  on  a  par¬ 
ticular  component  is  indicated  by  a  dash.  Events  with  a  SNNR  less  than  6.0  on 
the  bandpassed  data  were  not  used  in  the  computation  of  the  average  improve¬ 
ment  and  are  indicated  by  an  asterisk.  Chirp  lengths  which  are  not  considered 
to  be  optimum  are  also  indicated  by  an  asterisk  and  likewise  are  not  included 
in  the  computations. 

The  majority  of  regions  have  an  average  from  1  to  Z  dB  SNNR 
improvement  with  matched  filtering.  The  exceptions  are  the  Eastern  Kazakh 
test  area  which  have  poor  chirp  filter  improvements  (less  than  1  dB)  and  the 
Western  Pacific  area  from  Taiwan  to  the  Kuriles,  which  have  improvements 
above  3  dB.  The  results  of  these  data  differ  little  from  Special  Report  No.  7 
except  in  the  Taiwan  region  where  re -analysis  of  several  previous  events  has 
improved  the  average  gain  in  that  region. 

Omitting  those  events  with  a  bandpassed  SNNR  of  less  than  6, 

87  percent  of  the  Love  wave  chirps  have  positive  SNNR  improvement  and  45  per¬ 
cent  have  Z  dB  or  greater  improvement.  For  the  Rayleigh  wave  91  percent  of 
the  chirps  have  positive  SNNR  improvement  and  50  percent  have  improvements 
greater  than  Z  dB.  At  ALP  A,  60  percent  of  the  Rayleigh  wave  chirp  improve¬ 
ments  exceeded  Z  dB.  A  Z  dB  increase  in  SNNR  is  equivalent  to  about  an  in¬ 
crease  of  0.  1  m,  units.  For  the  1 97Z  winter  events  (Section  VI),  for  example, 

b 

this  increase  in  apparent  magnitude  would  raise  events  at  the  50  percent  de¬ 
tection  level  to  the  6Z  percent  level  and  events  at  the  90  percent  detection  level 
to  the  93  percent  level. 
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575.4 

335.1 

270.0 

2.3 

8.7 

3.6 

TA8LE  V-i 
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RcGION  *  T!  JR 

•  . 

GRFFC E 

, turkfy 

.BLACK 

SEA.F  FUPPPF t F 

Mcnj  JPRRAN 

c  AN 

EVENT 

CHIRR  LENGTH 

8AN0°A  5SED 

SNNR 

SNNR 

IMP 

DPI 

LOT 

LRV 

IRR 

LOT 

l  RV 

LRR 

LOT 

1  RV 

1.  CP 

TIJR  +  206+in 

100 

320 

320 

34.? 

34.4 

1°.? 

l.o 

3  .6 

7.3 

CRF*235*02 

?50 

?00 

200 

7.3 

8.8 

5.9* 

0.3 

3.2 

?  .6 

CrF*24?*22 

260 

?30 

?30 

11.4 

12.4 

R.2 

1.5 

3.7 

3.7 

C.RE*310*19 

?00 

40C 

400 

34.2 

53.6 

2°  .  ° 

1  .8 

6.5 

6.7 

C.PE*311*0O 

150 

?60 

- 

6.4* 

3.2* 

— 

0.4 

— 

CRE*321*03 

— 

350 

— 

- 

3.8* 

— 

— 

2.0 

- 

CRS*244*14 

160 

- 

- 

342.8 

- 

- 

3.0 

- 

- 

CYP*215*15 

150 

300 

?oo 

00.3 

14.6 

c 

• 

-t 

r\ 

1  .9 

1.5 

0.0 

onn*2lo*io 

180 

330 

- 

1R.5 

37.3 

- 

2.8 

-.6 

- 

GPE+2‘*1*08 

100 

3?0 

3?0 

36.0 

?5.6 

30.6 

1  .4 

4.0 

5.0 

GRE*23?*u4 

140 

300 

200 

8.8 

1  2.? 

4.3* 

5.0 

?.P 

GR  E*3  25*0° 

?60 

300 

200 

24.0 

1  8.0 

20.1 

0.0 

4.1 

4.3 

GRF*329*0l 

mo 

?60 

?  60 

52.9 

35.9 

30.  n 

3.1 

3.1 

3.1 

GR  E*329*03 

180 

?6C 

260 

145.7 

394.3 

144.0 

4.0 

2.3 

3.1 

PRF*330+15 

220 

260 

?  30 

4.  8* 

4.1* 

5.8* 

0.7 

4.1 

1.7 

MEO*219*03 

150 

?  50 

250 

2?. 4 

15.3 

15.0 

-1.1 

4.1 

7m 

MFD+242*2 1 

200 

no 

350 

31.7 

37.8 

35.3 

1.0 

3.8 

4.0 

MFD*320*12 

150 

350 

10.1 

3.6 

7.3 

3.6 

1  .o 

0.4 

MFO*333*13 

350 

400 

10.6 

26.6 

20.8 

5.2 

0.7 

0.3 

M  F0*334*0 1 

- 

?00 

— 

3.4* 

— 

— 

-.3 

— 

TUR  *2 16*0? 

’00 

200 

200 

1  92.3 

67.7 

49.4 

1  .2 

0.7 

1  .9 

TlJR  *216*21 

130 

200 

200 

309.1 

151.3 

10?.? 

1.5 

1  .0 

-.! 

TIJR*2!7*05 

160 

?40 

- 

1  26.6 

*3.1 

- 

1.1 

1  .9 

- 

TIJR*220*05 

130 

3?0 

240 

14.7 

14.4 

3.0* 

1.5 

n  .4 

3.1 

TIJR*236*21 

130 

240 

280 

10. P 

5.?* 

5.3* 

0.8 

4.9 

4.4 

WRS*?04*05 

80 

120 

— 

33.4 

58.4 

— 

1  .4 

0.5 

- 

=■  l  .on  7.59  ?.*■> 
=  1.36  1.01  l.°4 

=  5-*  44  77 


MOAN  SNNR  IMPROVEMENT 
STANP6Rn  PFVIATTPN 
MHMRFo  pr  FVFNTS 


TARIE  v-l 


CHIRP  MATCHED  FILTER  IMPROVEMENT  T  N  SIONAL+NDI  SF  /NE>1  SE  PATIO 

(DADE  4  np  ift) 


PVT, ION  ‘CAS 

•  — 

CASPIAN 

SEA, 

E  A  S  T  F  0  At 

c  alt  Asus 

EVENT 

CHIRP  l  ENDTH 

RA  NDn  A  S  SFD 

SNNP 

SNNP 

]  MP 

(DP) 

t  OT 

L  P  V 

|  PR 

l  0T 

l  RV 

|  RP 

l  0T 

1  PV 

l  oo 

C  AS / t  3  5/04 

?R7 

6° 

to 

0.0 

0.0 

0.0 

>.4 

1.2 

1  .0 

CAU/262/06 

- 

302 

- 

- 

0.0* 

— 

— 

0.0 

— 

C.AU/2RR/17 

30  F 

45* 

45* 

0.0 

0.0 

0.0 

2.F 

’.3 

"*  .4 

WK7/356/0 6 

241 

1  74 

I  76 

n.o 

O.o 

0.0 

1  .o 

1.4 

1.2 

C  At|*166  +  on 

7*0 

•»on 

300 

22.6 

74.4 

44.0 

0.7 

1  .3 

1  .0 

I OA*01fl*?1 

1?0 

600 

600 

43.1 

1  6.  7 

13. ■» 

3.0 

1  .  F 

MEAN  SNNP  TMPPnVFMPNT  *  ’.FT  7.03  1  .67 

STANDARD  DEVIATION  =  l.F?  1.13  1.D’ 

NIIMPEp  np  EVENTS  *  F  F  * 


V-IO 


TfPir  v-1 


I 

I 

I 

I 

I 

I 

I 

I 


CHIRP  M  AT  f  urn  FU.TFR  !  MPR  H  Vp  W  NT  *N  RlT,MAl+NnTcF/MPTSC  DAT!0 

(OAT,F  f  nc  IF) 


RFGTHM  'IRA*  -  FOIITHFCN  IRAN 


FVTMT 

third  length 

n  AMO  pa  c  <;p  n 

S  MM  R 

RMNP 

IMP 

(OR) 

1  ot 

L  R  V 

LPR 

1  OT 

1  RV 

|  RP 

LOT 

I  op 

L  DO 

! o  A*006*0° 

7  00+ 

420 

_ 

47.4 

71.1 

- 

-7.4 

0.0 

- 

IRA*014*?7 

7  50* 

560 

560 

n1,A 

77.4 

1  3.P 

-7.1 

-.0 

7.1 

!  R  A*041 +09 

700 

440 

440 

71 

7.7 

6.6 

-.0 

1  .7 

7.6 

I&  A*OMM  6 

A  50 

400 

400 

1  1  .A 

77.6 

«.K 

0.0 

o.o 

1  .0 

TR A*06B*?1 

750 

- 

400 

1  °  .  4 

— 

72.o 

0.7 

- 

7.7 

TR  A*1 56*0fl 

150* 

760 

- 

11.3 

7.  0* 

- 

0.4 

1  .7 

- 

TR  A*1  66*07 

- 

540 

*40 

- 

7.3* 

4.7* 

- 

5.6 

0  .6 

IP A*1S7*11 

1  60* 

700 

7  60 

46.4 

R.p 

«.R 

-7.7 

1.7 

-1  .1 

TR  A+167MO 

160* 

4R0 

4  PO 

63.6 

70.7 

16.0 

-1  .6 

7.7 

1  c 

• 

I P  A  *  1 64*  T 3 

100 

750 

+  60 

767.  0 

6  R.  4 

74.7 

0.6 

2.7 

1  .6 

IP  A*1  66*00 

160 

640 

4  p  0 

444.  0 

14R.0 

1  70.0 

-.4 

7.4 

7.C 

IPA*I  6R*?3 

160* 

700 

700 

7  7 . 7 

74.4 

0.0 

-7.4 

0.7 

1.7 

TR  A*1  76*00 

160* 

700 

+  00 

100.4 

?6.7 

71.3 

-.0 

7.0 

7.0 

TR A*1P4*12 

700+ 

- 

- 

1 777.7 

- 

- 

-1  .6 

- 

- 

TR A*1P4*14 

00 

760 

300 

?.n* 

7.4* 

7.6* 

1  .7 

7.7 

0.7 

I R  A  * 1 R5*07 

_ 

700 

700 

_ 

77.6 

71.6 

— 

7.7 

3  .4 

T  R  A*  ? 

- 

+  60 

- 

- 

7.7 

- 

- 

0.6 

- 

T  R  A*1 R5*7 1 

100 

760 

700 

747.6 

or 

• 

r- 

cr 

°3.7 

1  .3 

4.2 

■j  n 

• 

IR  A  +  l  ft  7*16 

T  00 

7  50 

?ro 

16.6 

21  .R 

71  .4 

3.0 

2.4 

1  .’ 

IR 6*1*7*21 

*  70 

- 

- 

7.7* 

- 

— 

1  .4 

- 

- 

IR  A*1  RP*06 

— 

+  50 

3  c  0 

- 

6.1 

6.7* 

— 

1  .6 

0.7 

T  R  A  *  1 03*  ?  + 

- 

700 

700 

- 

64.7 

67.0 

- 

7.7 

7.0 

I  R  A  + 1  06*  I  ■> 

1  on 

400 

- 

’1.7 

7.3 

- 

0.0 

0.*- 

- 

T R  A*  1  06*1  7 

100 

- 

- 

6.1 

- 

- 

1 .1 

- 

- 

IRO*060*0° 

7*n 

400* 

400* 

1  l.7 

3.6* 

4.4* 

1 .1 

1  .7 

-7.4 

I R  A*  ?  1 6*  22 

1  00 

700 

- 

4.7  + 

11.6 

- 

7.0 

-.o 

- 

I  R  A*21 7*0n 

+  0 

700 

7  00 

6.6 

1  7.7 

11  .c 

0.4 

o.l 

0.7 

TRA*?10*07 

70 

370 

7  PO 

c6. 6 

1  1  0.7 

67 . 0 

1.7 

0.6 

7.0 

To  A*?21*1Q 

70 

760 

o  60 

1  71.4 

4  6.0 

77.4 

7.4 

-.1 

4.7 

T  o  a*  720*06 

160 

700 

7  CO 

1  7.0 

17.4 

0.6 

7.1 

0.4 

0.4 

IR  A*3?l *10 

inn 

+  «0 

4  40 

14.4 

R.7 

°.7 

7.6 

1  .6 

7.0 

I RA*7?5*01 

160 

440 

+  on 

1  0.4 

4.6* 

S.o* 

3.7 

7.1 

7.0 

I R  A  *  3  30*2  ? 

170 

3P0 

7  PO 

60.  R 

10.4 

7  4.4 

3.1 

2.7 

1  1 
• 

MFAN  SMMP  I MPDOVr  Mr  MT  = 

1  .47 

1  .40 

7.0° 

5  T  A  VO  AR  0 

OF  VI  « TT  ON  = 

’.’I 

1.47 

1.47 

KJMMBP  R 

OF  rVFNTP  = 

17 

74 

V-1 1 


TABLP  V-l 


rHTPP  MATCHFO  CRTCR  l  MPRPV/F  Mr  NT  IN  SI  ANAL+NPI  SF/NDI'F  RATIO 

(PAOf  *  OF  16) 


RFOIPN  *PAK*  -  PAKISTAN,  FASTFRN  TRAN 


PVFNT  CHIRP  IF NOTH 

LOT  IRV  I  PB 


BANPPASScO  SNNP 
lqt  LRV  lrr 


SNNP  IMP  (OB) 
IQT  (  RV  I  PR 


(R A*0?P*0P 
IR  A*103*15 
PAK*167*1 1 
P AK*16?*1 1 
PAK*17q*10 

PAK*195*18 
P  A**231*10 


200  360 

220 

240  160 

TOO 

280  240 

240 

■*00  200 


•*60 

4.2* 

6.1 

3.9* 

-1.0 

4.2* 

- 

- 

2.9 

320 

26.1 

14. R 

14.0 

3.2 

260 

- 

4  6.6 

33.6 

160 

IT 

. 

or 

< 

132.1 

83.5 

0.7 

6.6 

- 

- 

200 

14.6 

27.8 

22.3 

0.7 

MR  AN  SNNP 

I  MPRPVP  MF  NT  * 

1.84 

STANOAPO  OFVIATIPN  * 

1  .45 

NIJMRT R  PF  FVFNTS  * 

3 

-1.2  0.9 

1.1  1.1 

2.3  1.4 

2.1  1.7 

2.8 

0.3  0.4 

1.73  1.16 
l.Kl  0.66 
6  4 


V  - 12 


TARLF  V-l 


CHIRP  MATTHm  PILTFR  IMpRO  VFMF  NT  IN  SI  GNA  L*NHT  SF /N01  SF  RATIO 

(PACF  7  OR  14) 

RFC, TON  *  T07  *  -  TAO/HIK,  K  I R  G I  7  f  HIN01J  KUSH,  N  AFGANISTAN 

fVFNT  CHIRR  LENGTH  RANPiPASSFO  SNMP  SNNR  IMP  (OR) 

|  QT  IRV  l  RR  LOT  I R V  LRR  LOT  I  °V  I  P0 

KRG/301/13  12?  3R7  381  0.0  0.0  0.0  l  .■>  ?  •°  l.p 

T0Z/147/00  348  30?  -  0.0  0.0  -  2.1  3.6 

T07/274/14  ?4fl  470  470  0.0  0.0  0.0  -3.2  3. 1  ''•q 

AFG*040*  1  9  ">50  -  400  4.4  -  13.4  3.8  "  U  • 7 

AFCi  *101+03  100  400  400  44.1  44.4  40. ">  0.8  -1.0  0.^ 

HNK*043*08  240  ‘>50  2*0  7.8  4.4* 

HNK  *  1 77*0?  TRO  -  -  48.7 

KRG*006*04  80  100  130  4.1*  °.  7 

TAO*077*00  400  340  340  40.3  1?4.3 

T  07  *005* l ?  440  450  440  4.7  4. 4 

HNK* 178*20  -  700  700  -  4.4* 

HNK* 170* 1 4  -  200  200  -  107.7 

AFG*2 1 5* l 2  110  340  -  2.8*  2.1* 

AFC,*??6*09  310  400  400  4.4*  0.7 

AFG*??0* 1 0  160  -  -  4.1* 

AFG*30R*23  100  440  440  24?.  Q  210.7  3.3  4.R  S.7 

AFG*320* 1 4  310  -  -  4.4  -  -  -1 - 

HNK*?23*0l  160  400  400  14.4  4.7*  2.1*  0.7  -.5 

HNK*?40*16  110  300  300  2.0*  ?.4*  ">.<)*  -.4  1  .«  ?.7 

T07*311*12  100  1  20  1  20  10.’  11.3  11.1  3.0  1.?  1  .Q 

MF  AN  SNNR  TMPRDVFMFNT  =  1.12  2. ’4  ’>.50 

OTANOARO  0FVIAT1HN  =  2.0°  1.01  1.41 

NUMRFR  CF  FVFNTS  =  1?  11  1 1 


3.7* 

2.1 

0.0 

7.0 

- 

-1  .4 

- 

- 

4.7 

-1  .* 

0.3 

2.2 

121.0 

1  .7 

1  .0 

l.7 

6.6 

2.3 

0.^ 

0.6 

4.0* 

— 

3.2 

2.4 

76.? 

- 

2.1 

2.4 

- 

1.2 

-.4 

- 

6.4 

0.4 

7.1 

2.R 

- 

?.4 

- 

- 

TAO|F  V-l 

fHTRO  M  A  T  C.  H  P  O  cfLTFQ  jMPpnvcHrNT  fM  GNM  +  Mn!  sr /MO!  <;r  oAT)n 

(PAr,F  p  OF  1ft) 

R F 0  ION  'TTR*  -  5  MNKTANG*  TJPCT,  N  PUP^A,  HMIAYAN  MOUNTAIN* 


TVFNT 

PHIRP  t  PNG th 

pa  Kino  a  F$Fn 

5NNP 

SNNP 

I  MP 

(OP) 

1  QT 

L°  V 

ipo 

LOT 

L&V 

!  OP 

LOT 

1  OV 

LPP 

S IN/ 21 0/ l 4 

174 

705 

305 

c 

• 

c 

0.0 

0.0 

’.1 

7.0 

0.7 

5IN/241  /I  r- 

- 

1  74 

- 

- 

0.  0* 

- 

- 

0.7 

- 

T  TR/123/00 

lPft 

‘  04 

7  05 

0.0 

0.  0 

o.n 

4.5 

1.’ 

7.4 

T  TB/155/70 

A  7  P 

47P 

0.0 

n.O 

0.0 

7.’ 

1  .0 

1  .4 

7  IR/302/1  7 

217 

R7 

07 

0.0 

0.0 

0.0 

7.7 

1  .0 

1  .’ 

PI|R*140*14 

ftOO 

700 

IP.  7 

’6.5 

13.7 

7.2 

-1.4 

-.4 

rm*iF4*ift 

■’•in 

4*0 

420 

4  2  •  0 

1  4.5 

I’.  R 

7.4 

7.P 

3.1 

1ND*1  54*?Q 

250 

440 

- 

10.« 

7.5 

- 

1  .o 

-.7 

- 

T 1R*075*06 

7on 

700 

’00 

40.7 

77.1 

74.0 

O.T 

7.0 

7.0 

TiP^iftn*77 

1  ’0 

1  00 

1  75 

1ft. 4 

5.P* 

4.7* 

1  .p 

1  .0 

7.3 

T  T  P*170*04 

240 

’ftO 

760 

?c  .2 

74.2 

14.3 

-.° 

3.7 

1 .1 

T  im>io«;*04 

200 

- 

- 

ft. 7 

- 

- 

7.2 

- 

- 

Tjq*ioq*o? 

9Cf) 

- 

- 

ftfr.4 

- 

- 

7.4 

- 

- 

T !B*108*07 

’00 

- 

- 

10. ft 

- 

- 

7.0 

- 

- 

7Ib+204M4 

150 

200 

’00 

54?  ft. 7 

<”4.4 

511. R 

’.P 

1.7 

-1 .0 

T ! R*704*?l 

150 

— 

- 

27.9 

- 

- 

7.7 

- 

- 

7  I  B*  70  5*  2  7 

1  50 

- 

- 

R.  ft 

- 

7.5 

- 

- 

T!P*?0ft*l'‘ 

1  5  0* 

’00* 

200* 

20. p 

2  0.7 

15.7 

-.7 

1.’ 

O.R 

r  j  n*?4 ?*? p 

i  50 

35n 

74Q 

17.5 

1  O.R 

1  1  .7 

0.7 

1.0 

-1.7 

5 TK  *734* 1 4 

400 

400 

400 

<4- 

• 

a 

7  4.Q 

IP. 4 

0.4 

-1.5 

-.4 

S  fK  *?->4*\R 

400 

’00 

— 

4.  P* 

4.7* 

- 

-.1 

2.0 

- 

5  IK**!  1*10 

•>O0 

200 

200 

4ft. 1 

7  0.4 

20.4 

0.0 

7.7 

4.’ 

MFAN  FNNP  IMnonvFMrfgT  -  ?.IQ  1  .  l.?4 

FTAMOAPn  ntviATTHN  =  1.44  I.™  1.P4 

MIIMRFP  pc  PVFNTS  =  IP  1  3  1? 


.! 

7 


V  - 1 4 


TAP  ir  V-l 


rHIRP  MATTHPr)  PIITPP  !MOpr)VFMrNT  |M  S  I  ONA  If  MHI  ST /MOf  S  F  D  ATI  n 

(PA r,f-  o  OF  1M 


prrinM  »ch! 

9  . 

r  FNTPAL 

CHINA 

rVFNT 

CHIRP  MOTH 

LOT 

L  0  V 

t  RO 

IS? 

470 

4  7<3 

CM  I  /??0/  17 

lr'7 

1  Of 

i  «S 

CMt/?«5P/n7 

174 

4  71 

4  71 

rH|*oiA*n7 

1R0 

sso 

rso 

CHTM«o*?3 

PRO 

— 

rHl*?0?*lS 

?10 

too 

6  00 

MON* l S3*!  1 

200 

- 

— 

mhm*?44+’ 7 

1  SO 

- 

- 

Y>  IN*  0S7* 1 q 

700 

700 

<S00 

rH!*??S*0? 

If  0 

?oo 

CH!*243MS 

ISO 

POO 

?co 

PA  Nn°A  SSP  0 

SNMP 

F  f'MR 

}  Mp 

(  0P| 

LO7 

1  PV 

t  PP 

t  QT 

1  PV 

L  RR 

0.0 

0.0 

0.0 

7.7 

t  .0 

4.7 

o." 

0.0 

0.0 

1  .4 

1  ." 

1  .4 

0.  07 

0.  0* 

0.0* 

?.S 

7.7 

7  .1 

3  7 . 7 

17.7 

0.7 

2.7 

1 .0 

7  7 

'  • 

tS.o 

- 

— 

1  .K 

“ 

3«.  1 

SO.  f 

70.  O 

c.« 

7.S 

0.  o 

- 

- 

1  .7 

- 

S.  1 

- 

- 

? 

- 

“ 

4**  .  7 

77.? 

70.7 

0.* 

-.7 

-1 

?  7  .  p 

0.7 

- 

1  .s 

-.P 

“ 

124.4 

114.7 

12P.7 

4.1 

0.4 

0.0 

MP  A  Kj  RMNP  TMPPnvFMTA'T  =  7.70  1  .e**  1  .ot5 

FTANOAPO  nrvTATjnM  =  1  .  S?  l.«7  7.?4 

M»|MRFP  nr  FVFNTS  =  10  7  f> 


V-l  5 


TARLF  V-l 


CHjap  MftTCHPn  PTtTPR  TMPRflVFMFNT  IN  6!  GNAL+NPT  SF /NTT  SF  AATTD 

( °  AGF  10  OF  161 


REGION  »S!N'  -  SOUTHW*  STFP  N*  WFSTFRN,  NDRTHWCSTFRN  SINK!  AND 


FVFNT 

GHTRP  LENGTH 

RAMOPASSFD 

SNNR 

SNNR 

IMP 

ID*) 

1  QT 

LAV 

LAP 

LOT 

1  RV 

LA  A 

LOT 

LAV 

L  P  A 

S  TN/166/2? 

87 

174 

174 

0.  0 

0.0 

0.0 

1  .6 

1.2 

0  .0 

SiN/166/22 

08 

217 

21  7 

0.0 

0.0 

0.0 

1  .4 

—  .6 

-.7 

1  .° 

SIN/170/17 

108 

217 

217 

0.0 

0.0 

0.0 

0.7 

1.3 

S IN/221/01 

104 

78 

7B 

0.0 

0.0 

0.0 

l  .0 

0.6 

0.7 

rPS/2‘>6/16 

113 

23® 

— 

0.0 

0.0 

2  .6 

2.6 

S  T  N/ 277/ 1 ? 

104 

174 

174 

0.  0 

0.0 

0.0 

1  .7 

1.6 

0.4 

PAT*058*22 

100 

300 

300 

7.7 

1  2.2 

8.4 

1  .2 

2  .7 

7  .6 

FK7  *078  +  07 

— 

50 

50 

- 

4.7* 

3.3* 

*• 

0.4 

1.3 

KRG*028*70 

80 

370* 

32  0* 

0.7 

4.5* 

7.7* 

0.8 

—  1.8 

-1  .3 
1.0 

STN*oo2*lo 

80 

100 

100 

70.4 

76.8 

67.4 

7.5 

1.1 

STN*04?*05 

100* 

100 

100 

53.3 

70.6 

62.6 

0.0 

1.1 

0.4 

STN*047*27 

80 

150 

150 

17.1 

“*1.8 

2’.4 

0.8 

0.6 

1.0 

S I N*064*04 

1  00 

200 

200 

10.3 

44.6 

10.1 

-.2 

1  .* 

0.7 

1.2 

1  .5 

8  TM*1P4*06 

80 

17* 

- 

20.0 

*.8* 

8  fV*lR7*ol 

40 

175* 

175* 

256.0 

78.2 

53.0 

1.1 

2  .7 

S  IN*1 87*04 

140 

175 

175 

8.0 

1  2.7 

8.4 

4.2 

3.2 

1  .6 

2.5 

1  .1 

S  TN*  1  «»?^1 0 

160 

17* 

175 

6.0 

35.6 

26.6 

7.2 

S  TN*200*03 

100 

- 

- 

7.4 

*■ 

0.7 

5 T  N*  2  75* 1 6 

60 

75 

76 

12.6 

11.1 

5.3* 

0.4 

1.7 

2.4 

1 .1 

S IN*  243* 1 7 

120 

75 

78 

1.0* 

2.6* 

2.8* 

2.° 

1.2 

STN*307*1 2 

60 

78 

1  7* 

3.8* 

2.6* 

4.1  * 

1  .0 

0.7 

-.6 

S  TN*716*1 4 

60 

75 

1  *0 

4.** 

4.6* 

7.4* 

-.5 

-.2 

1  .6 

l  .77 
0.°R 
14 

2.7 

STN*325*05 

60 

1  50 

24.0  3.5* 

MFAN  SNNR  1MPP0VFMF  NT  = 
STANDARD  DF VT  ATT  ON  * 
NMNBPR  OF  FVFNTS  = 

1  .1 

1.61 

1.06 

18 

0.07 

0.80 

l7 

V-16 


TABLE  V-l 


CHIRP  MATCHED  FILTER  ImrrhvFMF  NT  IN  SI  GNAL*-NPI  SF/NM  SE  R  AT  I  n 

IP8AE  n  OF  16) 


RCCION  '  FK7 '  -  EASTERN  KAZAKH  TF  ST  AREA,  CENTRAL  KA7AKH 


FVFNT 

CHIRP  LENGTH 

BANOPASSFD 

SNNR 

SNNR 

1  MP 

(OR) 

I  OT 

LRV 

LPR 

LOT 

1  RV 

L  RR 

LOT 

LRV 

LPP 

FK7/145704 

113 

— 

- 

0.0* 

- 

- 

4.0 

FK7/IS7/04 

01 

200 

200 

0.0 

0.0 

0.0 

1.3 

1.8 

0.2 
•  A 

FK7/I70/04 

S3 

130 

1  30 

0.0 

0.0 

o.o 

o.s 

0.0 

1 

FKZ/282/06 

- 

i  30 

- 

- 

0.  0* 

• 

7.7 

«  £ 

FK7/?°4/06 

■jq 

122 

12? 

0.0 

0.0 

0.0 

0.0 

0  •  p 

-1  «n 

FK7 /333/06 

S2 

60 

60 

0.0 

0.0 

0.0 

1  .4 

O.s 

i  .i 

EK7/364/06 

6S 

R7 

p7 

0.0 

0.0 

0.0 

1  .■* 

1.1 

1  .7 

FK  Z  *07 0*04 

so 

100 

100 

10.0 

0.0 

6.6 

1.0 

1  .4 

1.7 

FK7*0PR*04 

130 

SO 

PO 

2.  P* 

4.0* 

3.’* 

4.2 

0.1 

O.o 

EK7  *2  20*03 

- 

60 

PO 

6.8 

^  ♦ 

0.2 

0  • 1 

FK7*307*01 

40 

PO 

100 

21.0 

27.8 

20. S 

-.1 

l.P 

1.6 

FK7  *34S*04 

SO 

SO 

7S 

77.1 

•>4.  0 

70.1 

-.3 

o.o 

1  . 1 

WK7*233*Q7 

60 

OO 

7S 

11.2 

27.8 

77.  S 

0.1 

O.o 

1  . 1 

WRS*204*0S 

- 

- 

120 

- 

4 . 0 

" 

'  .7 

117  P*306*04 

- 

80 

120 

— 

S.  7* 

4.4* 

o.o 

0.0 

K  A7  *  18 1*00 

_ 

— 

ion 

- 

- 

4.7* 

- 

«  « 

1  .4 

CK  7-240-03 

- 

80 

- 

— 

1  3.6 

1.1 

EK  7 -047-05 
FK7-204-0 1 

40 

40 

PO 

- 

R.7 

S.  «♦ 

OS. 2 

- 

1.1 

-.4 

1.2 

- 

mc  a  N  SNNP  T  MPROVF  mE  NT  = 

0.77 

1  .01 

o.oo 

STANDAPO 

deviation  = 

0.62 

0.S7 

0.07 

NUMBER 

PF  FVFNTS  = 

10 

12 

10 

V-17 


TA8LF  V- 1 

fHTPP  maTTHFO  c  It  TPP  lMpRnvrHfF'T  IN  SI  ONAI  +  N01  SR /NO!  SF  PATIO 

(  P  A  fi F  1’  HP  lft  | 


RFC!  ON  *  T  A  I  '  -  TATWAM,  <5  RYUKYU  !S|ANOS»  FAST  CHINA  SFA 


FVPNT 

rwfpp  1 PNOTH 

f)4  NOPA  SSC  0 

SNNR 

SNMP 

1  MO 

(OR) 

1  OT 

1  PV 

OR 

t  OT 

LPV 

L  P  D 

LOT 

|  0\f 

LPR 

rHT*028*04 

230 

700 

7C0 

11.7 

7.5 

6.9 

1  .8 

0.0 

7.7 

r-  f\ 

RYU*1 55*0? 

TOO 

500 

500 

1  08.7 

101.7 

ftft  •  5 

1  .2 
7.1 

7  .ft 

5.0 

QYM*107*0? 

240 

570 

520 

24.7 

48.  0 

74.1 

7.0 

7.4 

1  .7 

PYU*209*lft 

760 

ft40 

440 

261.4 

1  Oft.  R 

68.5 

2  .9 

7.0 

1  .8 

T  A  I *00 A* OS 

500* 

1000 

— 

5.0* 

7.2* 

7.0 

T  A  1*004*! 2 

?80 

ftOO 

400 

97.0 

64.0 

50.7 

4.5 

1.5 
0.  1 

7.1 

T  AT  *006*06 

700 

ftOO 

600 

19.7 

10.5 

17.5 

7.1 

_.o 

T  A  T  *01 0*05 

700 

800 

800 

50.6 

70.0 

17.7 

0.2 

4.0 

7.8 

T  A  1  *  1 70*08 

700 

500 

600 

158.7 

04.  8 

81  .4 

2.7 

4.7 

4.7 

T A!*IR?*1P 

200 

500 

CC0 

548.0 

157. ^ 

PO.  1 

4.7 

0.8 

?  .4 

T  A  T* 195*  2  7 

TOO 

— 

- 

20.1 

- 

- 

4.0 

- 

4.1 

1.1 

TAT*19«*1 ^ 

700 

400 

400 

TO.O 

55.8 

7  5.4 

7.6 

2.1 

7.8 

RYU  *  2  7  1  *  1  P 

710 

770 

770 

27. R 

70.  8 

ft. 5 

1  .ft 

T  A  1 *235*2 1 

750 

- 

- 

7T.o 

— 

“ 

7.0 

C  rt 

TAl*312*0ft 

340 

440 

440 

!  Rl.P 

76.1 

74.0 

4.1 

5.5 

5  .O 

T  A 1  *  3  ?0* 08 

7  70 

?ftO 

260 

74.7 

75.8 

76.5 

7.7 

2.7 

1  .7 

TM*7'>ft*0? 

720 

440 

T  80 

1  R°. ft 

c8.1 

57  .  ° 

4.4 

5.6 

5.4 

HnN*70ft*04 

740 

2ftO 

26C  14.4  70.4  14. 1 

Mcam  SNNP  TMPPnVcM5NT  = 
STANOAPO  OF V! ATI °M  = 
MIjMprp  nc  FVpNTS  = 

2.7 

7.0° 

1  .47 

1  7 

2.6 

7.01 

1  .O* 

15 

7.5 

7.11 

1  .84 
15 
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TARtF  V-l 


(-  H  1  P  P  MATOHCp  cIlTt:o  TMpppVPMf  MT  TM  5  1  ONA  l  +-MP1  /  NO I  pATTn 

(  p  a  r, f  1  ’  nr  1  f>  > 

RcpIOM  *K  UR  *  -  KIJOILF  lStAN^S  JAPAN  T(1  KAMCHATKA 


py/PNT 

("HlRp  lfmoth 

RANOOAC^rO  SNNR 

SNN'R 

IMP  ( 

np  ) 

LOT 

1  RV 

LOR 

1  OT 

1  PV 

LOP 

1  or 

|  0 \l 

1  PP 

KUR/175/21 

fl?  R 

R04 

P04 

0.0 

0.0 

0.0 

-.4 

4.0 

4.1 

A  7 

KIJP/146/01 

564 

76? 

74? 

0.0 

0.0 

0.0 

•  7 

4  .  q 

i 

*  • 

-T  1 

KUR/ 147/16 

470 

R7p 

P?P 

0.0 

0.0 

o.o 

4  ,o 
c  o 

'  •  1 
/.  Q 

KIJo/  15?/?1 

56  6 

497 

4«7 

0.0 

0.  o 

0.0 

7.7 

/  t 

KMO/190/1.6 

q?R 

1002 

100? 

o.o 

0.0 

0.0 

4.0 

4  .*• 

KUR/191  /or 

741 

PSP 

95R 

0.0 

0.  0 

0.0 

4.7 

4.0 

7.2 

A  C 

KUR / 1 91 / 09 

R?3 

071 

P71 

o.o* 

0.0* 

O.o* 

0  .R 

0  .P 

t  n 

K1IR/213/0? 

4RO 

QQR 

poi 

0.0 

0.0 

0.0 

4.7 

P.7 

K  AM*078*1 fl 

9«0 

- 

- 

3.6* 

— 

1  .  R 

<1  O 

KIJR*001*1  4 

700 

- 

— 

’.4* 

1.7 

KUR  *001  MR 

— 

700 

- 

- 

4.5* 

- 

- 

0  .0 

/  A 

KiJP*005*02 
KUR*00o"M 4 
KUR*02?*01 

700 

Q00 

700 

QOO 

°co 

1  1  RO 

4.7* 

4.7* 

7.5* 

7 .4* 

■>  .  o* 

4.4* 

7  .4 

4.1 

1  .0 

.2 

7.0 

0.0 

KUR*07R*?7 

9P0 

— 

-*.6* 

-4  .  * 

KUR *04  4* 1 6 

R40 

R  SO 

PRO 

’.4* 

4.  R* 

?.°* 

0.5 

R.O 

_  “> 

t  *T 

KUR  *049*\ R 

°40 

RRO 

RRO 

4.7* 

4.o* 

7  .  P  * 

7  .R 

~.  R 

a  7 

1  •  f 

A  *1 

KUR*054*0’ 

420 

1000 

1000 

’4.  P 

?  4.  4 

7’.’ 

7  .o 

0.7 

A  A 

1  • 

KtJR*0S4*27 

T40 

1000 

1  ooo' 

75.0 

1  o.  4 

14.9 

7  .P 

7.  ? 

7  •  1 

A  O 

KMR*054*3  7 

400 

1000 

’  ooo 

’0.7 

74.4 

71.7 

’  •  1 

\  •  1 

J  •  * 

KIJR*055*  1 0 

980 

1000 

1000 

44.7 

174.4 

00.4 

7.5 

5.7 

1  c 

R.7 

A  *5 

KUR*056*2  7 

4  70 

700 

700 

’  .5* 

4.0* 

7.7* 

1.7 

-  1  .  s 

0  • 

KUP*057*0? 

QRO 

700 

700 

4  4.4 

1  3  7  .  S 

107.1 

7.7 

1  »p 

7  1 

n  #  r* 

-*  n 

KUR*057*05 

onp 

PRO 

PRO 

4.  P 

7.4 

4.4 

7  .4 

T  •  1 
y  1 

• 

1  Q 

KllR*067**>3 

740 

1000 

1  ooo 

4.  R* 

11.  R 

4.4 

’  .  0 

4  •  1 

»  '  1  9  ** 

KUR*070*04 

420 

RRO 

- 

7  .  R  * 

3.0* 

- 

7.7 

1  .R 

- 

KUR*077*07 

9R0 

- 

- 

11.4 

• 

4  .* 

A  A 

^♦1*57*00 

*500* 

700 

700 

R.R 

4.0* 

4.4 

7.6 

4.7 

7  #  ^ 

A  0 

KUR* 1 60*1 9 

740 

1000 

1  •’00 

U.7 

1  0.7 

R  .  7  * 

?  .5 

^  #7 

r)  •  M 

KUR  *171*1  0 

ORO 

700 

1  1  50 

14.7 

70.7 

°.  ? 

0.2 

0.4 

7  #o 

KUR  * l 7 l*  7  7 

1000 

- 

- 

q .  q  * 

- 

- 

’.7 
«  £ 

/  c 

KUR  * 190*  2 1 

ROO 

R  SO 

1  OOO 

7.7* 

R.O* 

’  .0* 

7 

I 

*7  n 

KUR  *1 9?  *  Of 

- 

«sn 

— 

• 

R1  .  1 

•  r » 

*7  /. 

KjjR+1  94*00 

R60 

Q  50 

1  ooo 

?  7 . 7 

57. R 

70.4 

4.5 

7.7 

7  7 

KUP*19S*1 r 

ROO 

70  0 

700 

4.  «* 

O.  Q 

4  .  R 

4.7 

**  .  7 

7  • 

V-19 


TARLF  V-l 


CHIPO  MATCHFO  F!LTFR  TMPROVFMFNT  1N  51GNAL+NOI SF/NPISF  °ATln 

(PAGF  14  OF  1#>) 


RFC  TON  •«»>*•  -  K I  IP  ILF  ISLANDS  PROM  JAPAN  TO  KAMCHATKA 


FVFNT 


KUR*21 1+21 
KMR*216*0? 
KUR+217*1 7 
KIJR*2?6* 1 1 
K1IR*229* 1° 

KIJR*?31*21 
KOO*?32*17 
Kl  I R +  232  +  21 
KUR+232+23 
KHR*?  33+03 

KIIR  +  236*  10 
KUR+306* 16 
KlJR  +  31 1*1  6 
KUR+3 1R+01 
KMR  +  ^2‘»*1S 

KM R *377*17 
KUP+323+17 


CHIRP  I.FNG  IH 


LOT 

LRV 

LRR 

500 

700 

700 

500 

R  50 

850 

7?0 

030 

ORO 

940 

9R0 

- 

330 

1120 

1120 

RftO 

- 

ftoo 

R40 

ORO 

ftl  0 

- 

- 

770 

QRO 

9P0 

- 

1170 

— 

720 

qflo 

980 

«00 

700 

700 

- 

700 

- 

610 

- 

— 

ftoo 

— 

• 

R30 

ORO 

- 

7?  0 

700 

700 

ranopassfo  SNNR 


LOT 

1.RV 

LOR 

LOT 

ftft.ft 

Rl.l 

34.6 

2.4 

10.0 

14.7 

11  .R 

l.R 

1  88.6 

70ft.  R 

1  2R.0 

7.0 

4.5* 

ft. 6 

- 

1  .0 

3.3* 

3.6* 

3.1* 

1.5 

1  O.ft 

- 

- 

1ft. 1 

13.3 

10.3 

ft. 6 

3.3* 

- 

- 

3.4 

R2  •  ft 

ft6  •  6 

46.8 

3. ft 

- 

7.9 

— 

• 

ft. 3* 

4.0* 

ft.l* 

-1.1 

21.4 

71.0 

0.0 

5. ft 

- 

11.2 

- 

— 

ft.l  ♦ 

- 

- 

0.4 

4.7* 

— 

• 

-.3 

3.1  * 

3.7* 

- 

0.9 

ft. 2* 

ft.  8* 

4.7* 

-.6 

SNNR  IMP  fDR) 
l  PV  L  PR 


1.7 

R.O 
7.0 
ft. 7 

7.7 

ft.  R 

6.6 

3.0 

1  .ft 

6.R 
ft. ft 
4.1 


O.R 

7.7 


2.P 
7.? 
ft. 3 

3.0 


1.7 

’.1 


1.1 

4.7 


ft.l 


MEAN  SNNP  I MPR OVF MF NT  *  ’.6ft  4.71  4.73 

STANDARD  PFVIAT1DN  =  1.91  7.33  7.6° 

nijmh cr  nF  FVRNTS  *  ’3  20  24 


TARLF  V-l 


CHIRP  MATCHFD  FILTFP  I  MPP  OVF  MF  NT  TM  ST  GNA l ♦NOT  SF/NOI SF  RATIO 

(  P A 0 F  1*5  °F  16) 


PFOION  *K  AM »  -  KAMCHATKA  °FMT  N  SHI  A  ,  KOMANOCPSKY  ISIANOS 


FVpNT 

CHIRP  1  FNC.TH 

LOT 

1  P.V 

t  PP 

KAM/166/14 

610 

1046 

1046 

K0M/14A/10 

523 

6Q7 

697 

K  AM*003*06 

560 

700 

700 

K  AM*007*19 

TOO 

- 

- 

K AM*004*07 

600 

900 

POO 

K  AM*005*  1  6 

620 

R50 

P50 

K  AM*009  *1 4 

620 

- 

- 

K  AM*012*20 

500 

10p0 

1200 

KAM*042*7i 

740 

100  o 

P50 

K AM*057*27 

— 

1000 

K  AM*156*07 

- 

1  1  50 

- 

K  AM* 1 57*04 

620 

1300 

700 

K  AM S*14 

- 

1000 

- 

KAM*177*17 

620 

1160 

- 

K  AM* 1 RO* 1 4 

— 

non 

1  000 

K  AM* 1 86*1 3 

- 

650 

- 

K  AM*1  89*0*5 

- 

1300 

- 

K  6M*1<?3*08 

500 

- 

- 

K  AM*lQ7*n 

- 

1000 

- 

KAM*706*17 

— 

700 

700 

K  0M*044*22 

74  0 

700 

700 

KOM*151*21 

QRO 

700 

700 

KUR*709*00 

*00 

1  150 

1  1  *0 

K  AM*716*03 

610 

700 

700 

K  AM*216*1 2 

610 

700 

700 

K  AM*?3i*19 

500 

P60 

P50 

K  AM*737*0R 

610 

R40 

700 

K  AM*  7 1 7*0 P 

500 

700 

°40 
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The  number  of  events  not  detected  on  any  bandpassed  beam 
but  detected  with  a  chirp  filter  was  15  (out  of  556)  at  ALPA  and  13  (out  of  515) 
at  NORSAR.  Ten  of  these  13  NORSAR  events  were  from  the  Kurile-Kamchatka 
area.  These  regions  tended  to  have  larger  chirp  gains  than  the  other  regions, 
therefore  this  preference  is  not  unexpected.  The  exclusion  of  these  events  from 
the  detection  threshold  estimation  would  not  change  significantly  the  90  percent 
detection  threshold  of  this  region  bit  would  appreciably  raise  the  50  percent 
detection  threshold  level.  This  is  discussed  further  in  Section  VI. 


SECTION  VI 


NORSAR  LONG-PERIOD  SURFACE  WAVE 
DETECTION  CAPABILITY 

The  detection  history  of  the  517  events  in  the  data  base  was 
used  to  estimate  the  NORSAR  long-period  surface  wave  incremental  detec¬ 
tion  probability  curve.  Smoothed  estimates  of  detection  probabilities  were 
obtained  by  a  new  method  (Ringdal,  1974).  This  method  models  detection 
as  a  random  Gaussian  process  and  obtains  maximum-likelihood  estimates 
of  the  process  parameters  based  on  the  experimental  data  at  hand.  For 
a  detailed  description  of  the  theory  and  its  limitations,  and  a  derivation 
of  the  likelihood  functions,  the  reader  is  referred  to  Ringdal  (1974). 

Detections  were  determined  independently  for  each  component 
except  for  the  1971  events  where  detections  were  assessed  only  on  the 
vertical  component. 

The  criteria  used  to  classify  the  events  as  detected  or  non- 
detected  were: 

•  A  peak  on  the  bandpass  or  matched  filtered  trace,  in  a  signal 
gate  defined  by  a  4.  0  km/ sec  arrival  and  a  2.  5  km/ sec  arrival, 
which  is  3  dB  larger  than  any  other  peak  within  a  1  5  to  30 
minute  gate  (depending  on  distance)  centered  at  the  peak. 

•  The  peak  and  its  code  should  show  signal-like,  i.  e.  ,  dis¬ 
persive  characteristics. 
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After  an  examination  of  a  number  of  events  from  a  particular 
region,  the  analyst  gains  an  intuition  for  the  output  trace  characteristics  that 
play  a  major  part  in  marginal  detection  cases.  In  most  of  the  marginal 
detection/non-detection  cases  the  judgment  for  detection  was  conservative  with 
the  result  that  the  false  alarm  rate  was  low  (estimated  at  less  than  1 
percent)  and  that  the  detection  threshold  was  also  conservative. 

The  subsets  of  events  used  for  detection  threshold  estimation 
and  their  thresholds,  based  on  the  assumptions  of  the  detection  model,  are: 

•  All  earthquakes  in  Eurasia,  Figure  VI-1.  The  90  percent 

detection  threshold  is  at  m  =  4.  5  which  is  the  same  as  reported 

b 

in  Special  Report  No.  7  for  NORSAR. 

•  Kurile-Kamchatka  area.  Figure  VI-2.  The  90  percent 

detection  threshold  is  at  rn  =  4.  5  which  is  less  than  Special 

b 

Report  No.  7  but  is  0.4  M  units  higher  than  ALPA.  ALPA 
is  substantially  closer  to  this  area,  however,  than  NORSAR. 

This  area  also  had  the  majority  of  chi rp -filter-only  detections. 
Although  the  90  percent  detection  threshold  is  essentially  un¬ 
affected  by  these  events,  the  omission  of  the  chirp  filter  detec¬ 
tions  would  raise  the  50  percent  detection  threshold  level  from 

m  =  3.  9  to  about  m,  =  4.  1. 
b  b 

•  Central  Asia  events,  Figure  VI-3.  The  90  percent  detection 

threshold  is  m  =4.6  w  lich  is  greater  than  last  reported, 
b 

These  set  of  events  is  different  than  the  last  report,  however, 
as  this  set  includes  only  events  of  short  chirp  lengths  (  <  250 
seconds)  of  the  Sinkiang,  Tadzhik,  Tibet,  and  China  areas. 
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FIGURE  VI-2 

NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS 
FOR  THE  KURILE-KAMCHATKA  ARFA 
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Because  of  the  number  of  points  falling  outside  of  the  90 

percent  confidence  interval,  this  estimation  is  not  as  reliable 

as  the  two  previous.  The  90  percent  detection  threshold  value 

at  ALPA  for  this  area  is  m  =  4.  4. 

b 

•  Presumed  explosions,  Figure  Vl-4.  The  method  of  detection 

threshold  estimation  used  here  is  subject  to  possible  large 

errors  when  the  number  of  events  is  small  as  in  this  set. 

Although  the  indicated  90  percent  detection  threshold  is  at 

m  =  5.  3  this  estimate  is  highly  unreliable, 
b 

•  Winter  and  summer  events,  Figures  VI- 5  and  VI-6  respectively. 
The  winter  90  percent  detection  threshold  is  at  =  4.5  and 
the  summer  90  percent  detection  threshold  is  at  -  4.  1 
These  sets  contain  only  the  1972  events.  The  summer  threshold 
estimate  is  not  as  reliable  as  the  winter  estimate  due  pri¬ 
marily  to  the  noise  level.  During  the  summer  when  the  noise 
level  is  low,  a  small  change  in  absolute  noise  level  produces 

a  large  change  in  the  equivalent  surface  wave  magnitude  of  the 

noise.  The  winter  detection  threshold  curve  in  contrast,  has 

a  steeper  slope  and  narrower  confidence  limits,  indicating  a 

stable  estimate  of  the  detection  threshold.  Section  IV  of  this 

report  shows  that  an  average  of  3  to  4  dB  improvement  using 

MCF  processing  can  be  expected  in  the  winter  at  NORSAR.  This 

would  be  a  change  of  0.  15  to  0.  2  M  units  in  the  detection  thres- 

s 

hold  at  NORSAR  in  the  winter. 

The  difference  of  0.  2  M  units  in  the  90  percent  detection  thres- 

s 

holds  for  Central  Asia  events  between  NORSAR  and  ALPA  seems  unusual  since 
NORSAR  is  30°  to  40°  closer  than  ALPA  to  this  region.  This  difference  can  be 
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explained  by  the  spectra  of  the  Central  Asian  signals  (Eyres,  1972).  The  sig¬ 
nal  spectrum  of  the  events  in  Central  Asia  at  NORSAR  tend  to  be  strongest  at 
the  high  end  of  the  passband  limit  of  0.  059  Hz  which  is  also  near  the  predom¬ 
inant  microseismic  peak.  Thus  by  using  the  passband  of  0.  025  -  0.  059  Hz, 
much  of  the  signal  energy  of  these  events  is  rejected  which  in  turn  may  raise 
their  detection  threshold.  The  spectra  of  these  events  at  ALP  A,  however,  do 
not  contain  as  much  higher  frequency  energy  and  thus  are  less  attenuated  by 
bandpass  filtering. 


SECTION  VII 


BEHAVIOR  OF  STANDARD  DISCRIMINANTS 


The  three  standard  surface  wave  discriminants  AL  versus 

AR  versus  m  and  M  versus  nri  were  measured  routinely  from  the  three  com- 
b  s  b 

ponent  bandpassed  beams  of  the  381  detected  events. 

A.  M  -  m,  MEASUREMENTS 

s  b 

The  surface  wave  magnitude  (M  )  was  computed  from  the 

s 

bandpassed  beams  by  the  following  formulas: 

M  =  log^-  +  1. 66  log  A  (  A  ^  25°) 

s  1 

M  =  log  A  +  log  A  +  0.  92  ( A  <  25  ) 

s  T 

where  A  =  maximum  peak-to-peak  amplitude  of  signal  in  milli¬ 
microns 

T  =  period  in  seconds  of  cycle  corresponding  to  A 
\  =  epicentral  distance  in  degrees 

Figure  VII- 1  shows  the  vertical  Rayleigh  (LR)  wave  M^  versus 

m  and  Figure  VII-2  shows  the  Love  (LQ)  wave  M  versus  m  .  In  both  figures 
b  so 

presumed  explosions  are  plotted  with  an  asterisk  and  earthquakes  with  a  circle. 

Complete  separation  between  classes  for  this  data  set  is  obtained 

by  the  M  -m  discriminant  for  both  Rayleigh  and  Love  waves,  the  minimum 
1  s  b 
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separation  between  earthquakes  and  presumed  explosions  having  the  same  m^ 

is  0.  3  M  units  for  the  LR  vertical  and  0.8M  units  for  the  LQ. 
s  s 

These  M  -  m,  results  differ  from  those  of  Special  Report  No. 
s  b 

7  because  not  only  have  an  additional  138  events  been  added,  but  Love  wave 

M  values  have  been  computed  for  the  1971  presumed  explosions.  In  addition, 
s 

the  eight  events  previously  having  poor  separation  have  been  reinvestigated 
and  found  to  contain  some  mistakes  in  the  Mg  calculations  and  some  incorrect 

m.  values, 
b 

The  presumed  explosion  labeled  '1'  in  Figure  VII-2  is  actually 

a  double  event  (EKZ*345*04)  and  has  a  higher  than  expected  Mg  value.  The 

"PDE"  bulletin  did  not  assign  a  magnitude  to  the  second  event  thus  we  have 

used  the  m  of  the  first  event.  Results  from  the  short  period  NORSAR  indi- 
b 

cate  that  the  second  event  is  actually  larger  and  is  near  mb  =  6.  0  (Ringdal, 
1973). 

The  earthquakes  lying  closest  to  the  presumed  explosion  popu¬ 
lation  were  investigated  for  regional  biases,  radiation  pattern  effects  or  com¬ 
putation  errors.  In  particular  those  events  processed  also  at  ALPA  (Special 
Report  No.  8)  were  used  to  corroborate  the  NORSAR  values. 

There  was  no  region  which  had  consistently  low  Mg  values.  The 

events  in  common  at  NORSAR  and  ALPA  with  low  Mg  values  had  low  values 

at  both  arrays,  except  where  a  radiation  pattern  effect  was  obviously  the  cause, 

for  example,  as  in  YUG*243*00  which  had  Mg  Rayleigh  =  3.  0  and  Mg  Love  = 

2.  9  at  ALPA  and  had  M  Rayleigh  =  2.0  and  M  Love  =  2.  7  at  NORSAR. 

s  s 

Four  events  with  m,  values  taken  from  the  PDE  bulletin  caused 

b 

some  problems  due  to  averages  of  a  wide  range  of  mb  values.  For  example, 
YUG*180*01  had  an  average  m^=  4.  6  from  four  stations  with  >  15  .  The 
mb  values  were  4.  0,  4.  1,  4.  1,  and  6.  1  with  the  last  value  given  by  a  station 
over  80°  away.  Fortunately,  these  instances  are  rather  rare.  While  we  have 
uncritically  used  published  mb  values,  except  for  recomputing  where  possible 
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to  include  only  teleseismic  stations,  we  are  aware  that  careful  examination 
of  the  short  period  data  would  result  in  better  values  and  probably  in  better 
class  separation.  These  problems  are  discussed  in  Te::as  Instruments  Special 
Report  No.  9,  1973. 


AL-AR  VERSUS  m. 


The  AR  parameter  was  introduced  by  Brune,  Espinosa  and 
Oliver  (1963)  as  a  measure  of  the  total  Rayleigh  wave  energy.  The  AL  parameter 
is  the  equivalent  measure  of  the  Love  wave  energy.  Evernden  (1969)  has  used 
these  parameters  to  discriminate  between  earthquake  and  presumed  explosions. 
Our  AL-AR  measurements  were  made  from  the  bandpassed  beams  in  a  signal 
arrival  gate  corresponding  to  a  velocity  window  of  4.0  to  2.  5  km/ sec  and  have 
been  scaled  (Harley,  1971)  so  as  to  permit  a  comparison  with  Evernden's  results. 
The  scaled  AL-AR  values  were  then  normalized  to  a  body-wave  magnitude  of  5.0 
and  a  distance  of  20°  (Brune,  1963). 

Figures  VII-3  and  VII-4  show  AR  versus  and  AL  versus 

m  respectively.  In  both  figures  the  earthquakes  are  plotted  as  circles  and  the 
b 

presumed  explosions  as  asterisks.  (The  presumed  explosion  labeled  1  is  the 
double  event  on  December  10,  1972  (EKZ*34  5*04)). 

Both  AL  and  AR  isolate  the  two  populations  with  AL  being  slight¬ 
ly  superior  to  AR.  A  straight  line  at  AL  =  8.0  would,  with  the  exception  of 
EKZ*345*04,  give  complete  separation  of  populations.  A  line  at  AR  12,  the 
Highest  value  of  a  presumed  explosion,  would  include  two  earthquakes  in  the  pre¬ 
sumed  explosion  population.  Recomputing  AL  and  AR  for  EKZ*345*04  using  an 

m  =6.0  would  improve  separation  considerably,  however.  In  the  range  above 
b 

the  estimated  50  percent  detection  threshold  for  presumed  explosions  m  =  5.  1, 

b 

both  of  the  discriminants  are  good. 
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FIGURE  VII- } 

AR  -  m,  DISCRIMINANT 
h 

COMBINED  1971  AND  1972  EVENTS 


VII -6 


M  * 


°o  o 

o  0  ° 

0°  □□ 

0  00  80 

nnBDQ  Q  §0 


g°g§8 


O  o 


io  °0  °  °g0 

On  ^°8  ^00  0 

SS°°  §go0°  ° 


080  0fi  °ooe 

§“§I8  Yo  Jo  o  °S 

°  8Soo8b°0°oS§o° 

8rJb88i:§!os« 


O  8s  Jo§o° 
8  8  §8  § 
0qo®o  8° 


*  * 


In  conclusion  the  measurements  of  M  ,  AL  and  AR  versus  m 

s  b 

provide  good  discrimination  capabilities  at  NORSAR.  In  most  cases,  these 

discriminants  would  perform  even  better  if  better  m  estimates  were  available. 

b 

Careful  selection  of  available  information  with  greater  weight  being  given  to 
the  LASA  and  NORSAR  values  would  be  a  direct  approach  to  this  problem. 
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SECTION  Vlll 


CONCLUSIONS 


A.  MAJOR  RESULTS 

The  major  conclusions  from  the  five  studies  of  the  NORSAR 
Long-Period  Array  Evaluation  are  summarized  below.  These  conclusions 
are  based  on  all  the  results  obtained  during  the  past  two  and  one-half  years. 

1.  Noise  Analysis 

•  The  noise  level  in  the  20-40  second  band  is  strongly  seasonal. 
Summer  levels  (March  to  October)  are  around  7  m/r  RMS  and 
fairly  constant.  Winter  noise  is  more  erratic  with  levels  vary¬ 
ing  from  8  to  IS  m/i  RMS.  The  maximum  observed  level  was 
79  m/i  RMS. 

•  The  high  winter  noise  levels  are  carsed  by  storms  in  the  North 
Atlantic  Ocean  and  many  last  from  12  to  36  hours.  The  noise 
generating  mechanism  is  not  understood. 

•  Noise  directions  are  also  seasonally  dependent  at  NORSAR  but 

are  not  strongly  correlated  with  noise  level.  Winter  noise  tends 

o  o 

to  be  strongly  directional  with  preferred  azimuths  of  0  ,  250  , 
and  280°.  Winter  noise  coherence  is  high  at  the  frequencies  of 
peak  power.  Summer  noise  is  generally  isotropic  often  with 
weak  easterly  direction  and  has  low  coherence  particularly  in 
the  signal  processing  band  of  17  to  40  seconds. 
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2. 


Signal  .Analysis 

Based  on  energy  spectra,  uncorrected  for  instrument  response, 
the  majority  of  the  signal  energy  generally  lies  between  periods 
of  17  to  40  seconds.  Some  Central  Asian  events,  however,  con¬ 
tain  significant  energy  down  to  13  seconds. 

The  arrival  azimuth  of  the  signal  is  along  the  great  circle  path 
with  only  small  deviations.  Some  events  have  shown  shifts  in 
arrival  azimuth  deep  in  the  coda  and  at  the  higher  frequencies. 

Signal  similarity  usually  is  very  high  along  the  propagation 
direction  but  degrades  rapidly  normal  to  the  propagation  path. 

3.  Array  Processing  Performance 

a)  Full  Array 

•  Multichannel  filter  improvements  in  array  gain  are  sea¬ 
sonally  dependent  for  off-design  noise  in  the  signal  band 
(0.  025-0.  059  Hz).  The  MCF  averaged  4.  0  dB  more  ar¬ 
ray  gain  than  the  beamsteer  in  winter  noise  conditions. 
This  implies  that  the  use  of  an  MCF  in  winter  could  lower 
the  detection  threshold  by  0.  2  M^  units.  Little  gain  im¬ 
provement  was  obtained  in  the  summer. 

•  Signal  amplitude  degradation  of  the  beamsteer  processor 
is  2.  0  dB  for  the  Rayleigh  wave  and  2.2  dB  for  the  Love 
wave.  The  MCF  signal  degradation  is  slightly  greater 
than  the  beamsteer  for  the  samples  used. 

b)  Reduced  Array 

•  The  high  signal  similarity  over  the  reduced  array  sign¬ 
ificantly  increased  the  performance  of  both  the  MCF  and 
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beamsteer  processors.  For  the  22  noise  samples  inves¬ 
tigated,  the  MCF  averaged  4  dB  more  noise  suppression 
than  the  beamsteer  processor. 

•  In  contrast  to  the  full  array,  the  MCF  improvement  over 
the  beamsteer  is  strongly  correlated  with  the  beamsteer 
noise  level.  The  MCF  improvement  over  the  beamsteer 
for  the  reduced  array  is  slightly  less  se.isonally  depen¬ 
dent  than  is  the  full  array. 

•  The  signal  amplitude  degradation  for  the  beamsteer  is 
1.4  dB  for  the  Rayleigh  wave  and  1.6  dB  for  the  Love 
wave.  These  values  are  0.6  dB  less  than  for  the  full 
array.  MCF  signal  loss  is  slightly  more  than  the  beam¬ 
steering  loss. 

c)  MCF  Processing  Using  One  Noise  Sample  for  MCF  Design 

•  Using  a  single  sample  of  summer  noise  for  MCF  design, 
array  gain  improvement  of  the  MCF,  when  applied  to 
events  up  to  15  days  from  that  sample,  showed  no  cor¬ 
relation  with  time  separation. 

4.  Matched  Filtering  Performance 

•  Considerable  effort  was  expended  to  determine  the  signal  pro¬ 
cessing  band.  The  17-40  second  period  band  (0.  025-0.  059  Hz) 
seemed  to  show  slightly  better  overall  performance  than  the 
other  bands  investigated  and  was  adopted  as  the  routine  signal 
processing  bind  for  NORSAR. 

•  Chirp  filter  lengths  were  successfully  used  to  regionalize  the 
events.  Chirp  filters  giving  optimum  £>NNR  improvements 
were  found  to  have  lengths  which  vary  in  a  regular,  consistent 
fashion  over  the  Eurasian  continent  allowing  contour  maps  of 
equal  chirp  length  to  be  made. 
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5. 


The  SNNR  improvements  of  the  chirp  filter  averaged  2  dB 
more  than  the  equivalent  bandpass  filter.  For  a  given  region, 
the  standard  deviations  of  the  gains  were  typically  0.  6  times 
the  mean  gain. 

Reference  waveform  filters  (RWF)  were  not  used  on  a  routine 
basis  at  NORSAR.  The  marginally  superior  performance  of 
the  RWF  did  not  seem  worthwhile  because  chirp  filters  offered 
more  stable  and  uniform  results,  were  less  sensitive  to  the 
presence  of  signal -like  noise,  and  parametric  values  were  simple 
to  generate  and  determine. 

Chirp  filters  produced  a  small  increase  in  the  number  of  de¬ 
tected  events  at  the  50  percent  detection  threshold  level  but 
did  not  materially  affect  the  90  percent  detection  threshold 
level. 


Detection  Threshold  Estimation 

The  90  percent  detection  threshold  (with  a  false  alarm  rate 
estimated  at  less  than  1  percent)  was  estimated  for  the  various 
event  populations  using  a  maximum  likelihood  method.  Under 
the  assumptions  of  this  method,  the  90  percent  detection  thres¬ 
holds  for  the  events  studied  here  are: 

All  Eurasian  earthquakes  m  =  4.  5 

b 

Central  Asia  m  =  4.  6 
b 

Presumed  explosion  m  =  5.  1  (this  value  is  based  on  a 

b 

few  events  and  may  be  unreliable) 


Kurile-Kamchatka  areas  m  =  4.  5 

b 

Summer  events  m  =  4.  3 
b 

Winter  events  m  =4.5 
b 
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•  The  difference  between  the  summer  and  winter  90  percent  detecion 
detection  threshold  is  attributed  to  the  increased  noise  level  at 
NORSAR  during  the  winter  months. 

6.  Behavior  of  Standard  Discriminants 

•  The  M  -m  criteria  achieved  complete  separation  between 

s  b 

event  classes  for  both  Rayleigh  and  Love  waves. 

•  AL  and  Ar  also  achieved  complete  separation  between  classes 

and  appear  to  give  slightly  better  separation  of  classes  than 

M  -m  at  NORSAR.  At  ALPA  the  AL,  AR  versus  m  criteria 
s  b  ° 

performed  more  poorly  than  the  M  -m^  criterion. 

B.  SUGGESTIONS  FOR  FUTURE  ANALYSIS  OF  NORSAR 

The  results  from  the  NORSAR  evaluation  conducted  during  the 
past  two -and -one -half  years  have  pointed  out  areas  which  should  be  studied 
in  any  future  analysis  of  NORSAR. 

•  The  correlation  of  ambient  noise  level  with  the  presence  of 
winter  storms  in  the  North  Atlantic  ocean  suggests  two  areas 
of  investigation.  First,  the  noise  generation  mechanism  should 
be  identified  as  either  a  coastal  surf  action,  an  interaction  of 
deef'  water  waves  with  continental  shelves,  or  an  open  ocean 
mechanism.  This  would  require  fairly  comprehensive  wave 
and  weather  information  from  both  Norway  and  England  and 
possibly  seismic  data  from  English  stations.  Second,  multi¬ 
channel  filters  should  be  used  on  wintertime  data.  The>e  may 

be  fixed  or  possibly  time -adaptive.  If  the  first  study  demonstrates 
a  structural  effect,  fixed  MCF's  may  be  useful  for  processing 
of  winter  data. 

•  More  events  from  central  Asia  should  be  analyzed.  Some  areas 
shaw  unusual  behavior  in  both  spectral  content  and  dispersion. 
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particularly  around  the  Caspian  Sea,  Tadzhik,  and  Kirgiz. 

These  events  tend  to  have  particularly  large  high-frequency 
energy  content,  relatively  wide  dispersion,  compared  to 
Sinkiang,  China,  and  more  northerly  areas,  and  occasionally 
reverse  dispersion.  It  may  be  desirable  to  use  wider  processing 
bandwidths  for  central  Asian  events. 

Chirp  matched  filters  should  be  evolved  to  more  sophisticated 
non-linear  models  while  maintaining  simplicity.  This  may 
realize  si  mificantly  higher  signal -to-noise  ratio  gains  for  the 
central  Asian  events  in  particular. 

Maps  of  chirp  filter  length  were  successful  in  indicating  seismic 
region  boundaries.  There  is  some  evidence  that  the  boundaries 
of  the  tectonic  systems  of  Asia  can  be  correlated  with  these 
maps.  It  may  be  possible  that  such  maps,  computed  for  several 
seismic  stations,  may  be  helpful  in  pointing  out  geologically 
interesting  areas. 
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APPENDIX  A 


Tabic  A-l  presents  the  entire  ensemble  of  events  processed 
during  the  evaluation  of  the  NORSAR  Long-Period  Array.  It  includes  59  events 
from  1971  and  456  events  from  1972.  The  latter  group  includes  133  new  events 
processed  since  release  of  Special  Report  No.  7.  Processing  of  the  1971  events 
was  discussed  in  Special  Report  No.  5.  A  number  of  events  previously  included 
in  the  original  1971  and  1972  ensembles  have  been  deleted.  These  are  discussed 
in  Appendix  B. 

The  events  are  ordered  chronologically  and  are  identified  by  a 
name  consisting  of  a  three  letter  area  rode,  the  julian  date,  and  the  hour  of 
origin.  The  year  is  indicated  by  the  symbols  separating  the  three  parts.  A  "/" 
indicates  1971,  a  indicates  1972,  and  a  "  -  "  indicates  1973  . 

Othor  information  includes  the  origin  date  and  time,  epicenter 
coordinates,  delta  and  azimuth  from  NORSAR,  m^,  Rayleigh  and  Love  surface 
wave  magnitudes,  normalized  AR  and  AL,  depth  (if  reported),  and  coded  com¬ 
ments.  The  latter  consist  of  a  letter  indicating  the  source  of  event  informa¬ 
tion,  an  "X'  if  the  event  is  a  presumed  explosion,  and  the  number  of  sites  used 
to  form  the  array  beam.  The  source  codes  are  "P"  for  the  Preliminary  Deter¬ 
mination  of  Epicenters  (PDE)  bulletin  of  the  National  Oceanic  and  Atmospheric 
Administration,  *'LM  for  the  Seismic  Data  Analysis  Oenter  bulletin  for  the  LASA 
array,  ,rN '  for  the  NORSAR  bulletin,  and  "I"  or  "J"  for  the  Massachusetts 

*  Note  that  these  symbols  are  the  same  as  used  for  the  ALPA  evaluation 
but  are  displaced  by  one  year.  That  is,  for  ALPA,  a  "  -  "  was  used  for  the 
1972  events,  etc. 


A-l 


Institute  of  Technology  -  Lincoln  Laboratory  bulletin  for  the  International 
Seismological  Month  (ISM)  which  was  from  20  February  to  19  March,  1972. 


The  events  obtained  from  PDE,  and  those  from  the  ISM  group 

which  were  also  listed  by  PDE  were  checked  for  m  determination.  Any  m 

b  b 

values  reported  by  stations  closer  than  15  degrees  to  the  event  were  deleted 
and  the  average  was  recalcula 
was  placed  next  to  the  revised  m 


and  the  average  m  was  recalculated.  If  a  change  in  m  resulted,  an  asterisk 

b  b 


Non-detection  of  a  event  on  the  Love  or  Rayleigh  components 

is  indicated  by  a  dash  in  the  M  column;  a  "C"  in  those  columns  indicates  de- 

s 

tection  by  a  chirp  matchH  filter  only.  Love  wave  magnitudes  of  the  1971 
events  were  computed  only  for  the  presumed  explosions. 
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APPENDIX  B 


Appendix  A  does  not  contain  many  of  the  events  previously  in¬ 
cluded  in  the  ensemble  lists  of  Special  Report  No.  5  and  No.  7.  In  Report  No.  5,  of 
152  events  listed,  only  59  were  processed  through  matched  filtering.  Several  others 
were  used  only  for  measurement  of  signal  loss  from  beamforming  and/or  M 

s 

m  calculations  on  the  unfiltered  data.  The  remainder  of  the  edited  events  were 
b 

not  processed.  Because  of  the  large  number  of  events  which  have  been  fully  pro¬ 
cessed  through  matched  filtering  in  this  contract,  only  those  59  similarly  pro¬ 
cessed  events  from  1971  have  been  included  for  consistency. 

The  single  event  deleted  from  Special  Report  No.  7  to  the  present 

report  was  ITA*035*04.  This  event  was  reported  by  seven  non-teleseismic  stations 

(A  <  15°)  but  only  one  station  assigned  an  m  (m  =  4.8).  Since  the  event  was 

b  b 

not  detected  at  either  NORSAR  or  LASA  the  value  4.  8  was  considered  erroneous 
and  the  event  was  dropped. 


APPENDIX  C 


This  appendix  contains  the  lists  of  the  93  noise  samples  edited 
at  NORSAR  from  1  May  1971  to  27  March  1973.  The  data  are  arranged  in  the 
three  tables  chronologically.  Included  are  the  start  time,  duration,  and  the 
number  of  good  sites  available. 

Table  C-l  contains  the  25  noise  samples  edited  in  1971. 

Table  C-2  contains  the  55  noise  samples  edited  in  1972.  Table  C-3  contains 
the  13  noise  samples  edited  in  1973. 


TABLE  C-l 

1971  NORSAR  NOISE  SAMPLES 


Day 

Start  Time 

hr:  min 

Length 
hr: min:  sec 

Number  of 

available  sites 

1  May 

18:00 

0.  39. 24 

16 

1 5  Mav 

10:00 

01.33.52 

8 

21  May 

05:00 

01.  16.  48 

14 

30  May 

13:21 

0.  57. 28 

14 

13  June 

10:00 

04.  37.  20 

14 

20  June 

12:15 

01.  16.  48 

14 

10  July 

23:00 

0.  51.  12 

11 

20  July 

00.  59.  44 

18 

30  July 

16:37 

01.  12.  32 

15 

8  August 

20:20 

04.  24.  32 

19 

19  August 

19:  50 

02.04.  44 

17 

29  August 

17:00 

02.42.08 

15 

7  September 

20:37 

01.  55.  12 

16 

1  3  September 

18:00 

06.02.  40 

14 

18  September 

0  5:20 

03.03.  28 

16 

29  September 

17:00 

01.  50. 58 

17 

8  October 

0  5:  20 

01. 21. 04 

16 

18  October 

08:50 

0.  59.  44 

16 

28  October 

00:06 

01.  38.08 

16 

6  November 

04:00 

02.  20.  48 

14 

17  November 

20:15 

04.09.28 

20 

28  November 

06:00 

04.03.  12 

19 

6  December 

07:45 

02.  16.  32 

18 

14  December 

04:00 

02.  16.  32 

15 

27  December 

20:00 

02.  29.  10 

17 

TABLE  C-2 


1972  NORSAR  NOISE  SAMPLES 
(PAGE  1  OF  2) 


Day 

Start  Time 

hr:  min 

Length 
hr: min:  sec 

Number  of 

available  sites 

5  January 

13:00 

6.  24.  00 

14 

5  January 

19:00 

6.  24.  00 

14 

9  January 

16:30 

6.  19.  44 

18 

1 3  January 

8:30 

3.  35.  00 

13 

19  January 

4:30 

4.  45.  52 

14 

30  January 

6:00 

5.  49.  52 

18 

10  February 

18:00 

4.  37. 30 

19 

18  February 

5:30 

0.  21. 08 

20 

26  February 

7:30 

6.  24.  00 

15 

8  March 

|  ■ 

2.  54.  56 

16 

28  March 

7.06.  40 

18 

8  April 

10:30 

5.02.  56 

19 

23  April 

6:00 

4.  45.  52 

22 

17  May 

5:00 

1.  25.  20 

16 

27  May 

15:00 

3.07.  44 

21 

3  June 

21:00 

4.  16.  00 

17 

1 5  June 

3.  35.  00 

21 

7  July 

2.  42.  08 

16 

13  July 

2.  50.  40 

20 

23  July 

2.  29.  30 

21 

1  August 

3.03.  28 

16 

7  August 

02:00 

06.  02.  40 

22 

1  5  August 

11:30 

05.02.  56 

16 

18  August 

20:00 

0.  42.  40 

19 

20  August 

17:00 

07.06.  40 

21 

24  August 

22:30 

05.  02.  56 

17 

28  August 

1:00 

0.  42.  40 

21 

31  August 

08:30 

05.02.  56 

18 

6  September 

10:00 

0.  42.  40 

18 

9  September 

07:00 

07.06.  40 

21 

1 2  September 

17:30 

06.  02.  40 

17 

C-3 


TABLE  C-2 


1972  NORSAR  NOISE  SAMPLES 
(PAGE  2  OF  2) 


Da 

Start  Time 

hr:  min 

Length 
hr:  min:  sec 

Numbei 

available 

14  September 

01:10 

04.  11.  44 

17 

17  September 

0.  42.  40 

19 

19  September 

05.02.  56 

20 

24  September 

09:00 

04.03.  12 

19 

28  September 

14:00 

0.  42.  40 

19 

1  October 

14:00 

03.  33.  20 

22 

6  October 

04.03.  12 

20 

12  October 

04.03.  12 

19 

17  October 

10:00 

0.42.  40 

19 

18  October 

04:00 

07.06.  40 

13 

27  October 

03:00 

04.  03.  12 

18 

27  October 

22:00 

0.42.  40 

19 

4  November 

12:00 

04.  33.04 

22 

6  November 

8:00 

0.  42.  40 

20 

12  November 

21:00 

06.  06.  56 

16 

15  November 

23:00 

0.  42.  40 

21 

20  November 

00:00 

03.  33.  12 

18 

26  November 

20:00 

03.  33.  12 

20 

4  December 

04:00 

01.42.  24 

17 

6  December 

17:00 

0.42.40 

21 

10  December 

03.  33.  12 

22 

18  December 

■ftflwlifr  tiittBli 

06.02.  40 

20 

25  December 

03.  16.  40 

21 

26  December 

15:00 

0.  42.  40 

21 

! 

| 

i 


I 

* 

I 

I 

I 

I 

I 

I 

I 


TABLE  C-3 

1973  NORSAR  NOISE  SAMPLES 


Day 

Start  Time 

hr:  min 

Length 
hr: min:  sec 

Number  of 
available  sites 

1  January 

08:00 

06.  33.  20 

1 

19 

8  January 

16:00 

05.02.  56 

22 

16  January 

12:00 

06.  02.  40 

16 

23  January 

19:30 

05.  02.  56 

17 

31  January 

07:00 

05.  31.  52 

20 

6  February 

07:00 

03.03.  28 

19 

12  February 

06:00 

06.  02.  40 

19 

20  February 

00:00 

05.02.  56 

22 

28  February 

17:00 

04.  33.  04 

15 

8  March 

02:00 

04.  11.44 

19 

14  March 

15:00 

04.  11.  44 

22 

21  March 

06:00 

07.06.40 

22 

27  March 

21:00 

04.03.02 

22 
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APPENDIX  D 

SOME  COMMENTS  AND  OBSERVATIONS  CONCERNING  THE 
EFFECTIVENESS  OF  THE  USE  OF  LINEAR  GROUP 
VELOCITY  CURVES  IN  CHIRP  FILTERING 


The  chirp  matched  filters  which  have  been  used  routinely  on 
NORSAR  and  ALPA  long-period  data  to  enhance  signal  detectability  are  designed 
to  have  a  linear  change  in  frequency  within  a  given  passband  over  a  given  time 
interval.  At  NORSAR,  the  passband  has  been  0.025  to  0.059  Hz  (17  to  40  seconds) 
with  durations  typically  from  50  to  1200  seconds.  The  "best"  chirp  for  an  event 
is  the  one  which  produces  the  maximum  increase  in  peak  signal  amplitude.  Us¬ 
ing  a  linear  chirp  as  a  matched  filter  implies  the  assumption  that  the  signals 
being  matched  have  a  group  velocity  which  changes  almost  linearly  with  fre¬ 
quency.  This  study  compares  the  actual  group  velocities  of  a  small  set  of  events 
with  the  linear  group  velocities  implied  by  the  best  chirps. 

The  events  were  six  Italian  events  which  occurred  on  February 
4,  5  and  6,1972  in  close  proximity  to  one  another  and  which  had  high  signal-to- 
noise  ratio  (  >  10).  Table  D-l  gives  some  epicenter  information  and  Figure 
D-l  is  a  mapof  their  locations. 

The  Love  waves  of  these  events  are  nearly  identical;  visually 
the  time  series  overlay  one  another  almost  exactly  throughout  the  signal  dura¬ 
tion.  However  the  Rayleigh  wave  phases  are  not  as  uniform  with  the  two  events 
ITA*035*09NL  and  ITA*036*07NL  having  a  phase  change  in  the  middle  of  the 
signal  which  are  not  evident  on  the  remaining  four.  The  events  were  beam- 
steered  to  the  expected  signal  arrival  azimuth  using  a  velocity  of  4.0  km/ sec 
for  Love  waves  and  3.  5  km/ sec  for  Rayleigh  waves.  The  bandpass  filtering  and 
chirp  lilter  processing  were  performed  on  these  beamsteered  data. 


D-l 


EVENT  PARAMETERS  AND  OPTIMUM  CHIRP  LENGTHS 


The  highest  gain,  linear  group  velocity  chirps  for  these  data  have 
been  documented  in  Section  V  of  this  report  and  in  Special  Report  No.  7  and  are 
listed  in  Table  D-l  for  both  LR  vertical  and  LQ  transverse  components. 

The  measurements  of  the  actual  group  velocity  curves  were  made 
using  the  technique  of  band-pass  filtering  described  Kanamori  and  Abe,  (1968). 

For  both  modes,  the  seismograms  were  filtered  by  five  zero-phase  filters,  each 
having  131  points  and  using  a  2  second  sampling  interval,  which  had  these  pass- 
bands: 

0.020  •  0.030  Hz  centered  at  40  second  period 

0.027  -  0.037  Hz  centered  at  30  second  period 

0.035  -  0.045  Hz  centered  at  25  second  period 

0.045  -  0.055  Hz  centered  at  20  second  period 

0.054  -  0.064  Hz  centered  at  17  second  period 

The  time  of  occurrence  of  the  peak  amplitude  of  the  filtered  out¬ 
put  was  taken  as  the  arrival  time  for  the  center  period  of  that  filter,  and  from 
the  distance  (A)  listed  in  Table  D-l,  the  group  velocity  for  the  corresponding 
event  and  period  was  computed.  These  results  are  plotted  by  event  in  Figures 
D-2  and  Figure  D-3  for  the  LR  vertical  and  LQ  transverse,  respectively. 

The  maximum  amplitudes  of  the  filtered  traces  were  norma¬ 
lized  to  the  25  second  period  amplitude  for  each  event  and  are  plotted  in  Figure 
D-4.  The  relative  spectral  amplitudes  shown  in  Figure  D-4  help  to  explain  the 
large  variation  in  LR  group  velocities  at  17  and  20  second  periods  shown  in  Figure 
D-2.  The  events  with  low  spectral  amplitudes  at  17  and  20  seconds  are  probably 
contaminated  by  noise  because  the  noise  energy  is  higher  at  those  periods  than  at 
30  to  40  seconds,  thus  a  large  variance  of  expected  arrival  times  would  be  ex¬ 
pected.  High  spectral  amplitudes  for  the  same  filter  had  peaks  within  one  cycle 
of  one  another  giving  accurate  time  estimates.  Figure  D-4  also  shows  the  two 
different  LR  wave  sets  of  this  suite  of  events.  IT  A*035#09NL  and  IT  A#  0  3  6*  07  N  L 
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FIGURE  D 


RAYLEIGH  WAVE  GROUP  VELOCITIES  FROM  ITALY 
MEASURED  AT  NORSAR 
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FIGURE  D-4 

SPECTRAL  RATIOS  NORMALIZED  TO  25  SECONDS 


have  relatively  more  20  second  energy  and  le6S  17  second  energy  than  the 
others.  LQ  relative  spectral  amplitudes  are  more  consistent  between  events 
than  the  LR  data.  By  omitting  ;he  low  spectral  amplitudes  at  17  and  20  sec¬ 
onds  on  the  LR  vertical,  and  taking  averages  of  the  remaining  at  each  period, 
smoothed  LR  and  LQ  group  velocity  curves  were  obtained  and  are  shown  in 
Figure  D-5. 

Figure  D-6  shows  the  two  group  velocity  curves  of  Figure  D-5 
plotted  versus  frequency.  Also  plotted  on  Figure  D-6  are  the  implied  group 
velocity  curves  of  the  chirp  filters  used  when  routinely  processing  the  event 
ITA*035*17NL.  The  point  of  intersection  on  the  linear  group  velocity  curves 
at  0.  040  Hz  is  the  zero-phase  frequency  of  the  chirp  which  defines  the  time  at 
which  the  chirp  filter  output  maxima  should  occur.  Listed  beside  each  curve 
is  the  length  of  the  chirp  in  seconds.  The  maximum  gain  output  was  for  lengths 
of  75  and  65  seconds  for  LQ  transverse  and  LR  vertical,  respectively. 

Three  things  are  immediately  apparent  from  Figure  D-6.  First, 
the  recorded  group  velocity  curves  are  not  linear  with  frequency  over  this  band. 
Secondly,  none  of  the  chirp  filters  used  in  the  processing  fit  the  entire  curve. 
Third,  the  velocity  of  the  chirp  maximum  at  0.  040  Hz  does  not  correspond  ex¬ 
actly  with  the  average  0.  040  Hz  velocity  measured  with  the  narrowband  filter 
set.  Because  the  real  group  velocity  is  not  linear  and  c,ince  the  chirps  used  in 
the  processing  are  chosen  on  the  basis  of  maximum  gain,  then  the  velocity  im¬ 
plied  by  the  chirp  should  tend  to  approximate  that  portion  of  the  real  group  ve¬ 
locity  curve  that  contains  the  most  energy.  For  the  event  used  in  Figure  D-6, 
the  highest  spectral  amplitudes  are  at  20  and  25  second  periods  for  the  Love 
wave  and  at  17  and  25  second  periods  for  the  Rayleigh  wave.  A  shift  of  the 
linear  velocity  curves  for  the  chirp  of  maximum  gain  (65  seconds  for  LQ  and 
75  seconds  for  LR)  to  intersect  the  actual  group  velocity  curves  gives  nearly 
exact  fits  to  the  actual  group  velocity  curves  at  the  two  points  of  highest  spec¬ 
tral  amplitude.  However,  over  the  range  of  chirp  lengths  used  on  this  event, 
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the  chirp  filter  maxima  vary  by  less  than  1  percent  for  the  Love  mode  and  less 
than  10  percent  for  the  Rayleigh  mode,  so  the  chirp  lengths  giving  maximum 
gain  are  not  sharply  defined  at  least  for  this  case. 

The  arrival  time  of  the  chirp  peak  is  also  important  ^or  detec¬ 
tion  purposes.  The  chirp  peak  should  arrive  within  the  same  interval  that  the 
unfiltered  signal  should  arrive  which  is  (1/2.5  (km/  sec)-  1/4.  0  (km/  sec))  x  A 
(km).  A  non-linear  fr 3quency-time  chirp  that  was  a  better  approximation  to 
the  dispersed  signal  would  reduce  the  uncertainty  in  the  arrival  time  of  the 
zero-phase  energy  and  hence  aid  detection. 

In  this  report  and  in  Special  Report  No.  7,  the  optimum  chirp 
lengths  are  contoured  for  the  Asian  continent  for  LR  and  LQ  waves.  Since  it 
has  been  shown  here  that  the  optimum  chirp  length  is  a  function  of  the  received 
signal  spectrum  as  well  as  the  actual  group  velocity,  then  some  of  the  questions 
concerning  these  plots  can  be  resolved.  For  example,  since  the  group  velocity 
curve  is  an  average  over  the  entire  travel  path  from  source  to  receiver,  the 
rapid  change  of  chirp  lengths  around  the  Hindu-Kush  area  may  be  due  to  sig¬ 
nal  spectra  received  at  NORSAR  rather  than  extreme  velocity  changes.  Be¬ 
cause  the  actual  group  velocity  curves  for  Asian  events  at  NORSAR  are  un¬ 
known,  the  influence  of  the  signal  spectra  on  these  chirp  lengths  also  is  un¬ 
known  at  present. 

The  concept  of  a  non-linear  chirp  filter  is  similar  to  that  of  a 
reference  waveform  filter  (RWF).  Maintaining  a  file  of  chirp  parameters  for 
various  regions  would  be  like  maintaining  a  suite  of  reference  events.  The  dif¬ 
ferences,  however,  are  that  the  chirp  parameters  would  be  based  on  averaged 
information  and  wou'd  require  only  minor  amounts  of  storage,  and  that  the  re¬ 
sulting  chirp  would  have  constant  amplitude,  precluding  the  noise  sensitivity 
of  RWF's  discussed  in  Special  Report  No.  7.  Thus  a  non-linear  chirp  could 
be  expected  to  produce  more  gain  than  a  linear  chirp,  less  gain  than  a  RWF, 
but  with  greater  stability  than  the  RWF. 
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However,  even  though  chirps  and  RWF's  aid  detection  by  in¬ 
creasing  signal-to-noise  ratio,  these  gains  are  primarily  at  the  50  percent 
detection  level  and  not  at  the  90  percent  detection  level. 


